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^ . Abstract 

(N ■ 

We evaluate the next-to-next-to-leading order corrections to the hard-scattering amplitude of the 

^ photon-to-pion transition form factor. Our approach is based on the predictive power of the 
P-i! 

conformal operator product expansion, which is valid for a vanishing /3-function in the so-called 

conformal scheme. The Wilson-coefficients appearing in the non-forward kinematics are then 

entirely determined from those of the polarized deep-inelastic scattering known to next-to-next- 
• »-h , 

^ . to-leading order accuracy. We propose different schemes to include explicitly also the conformal 
symmetry breaking term proportional to the /3-function, and discuss numerical predictions cal- 
culated in different kinematical regions. It is demonstrated that the photon-to-pion transition 
form factor can provide a fundamental testing ground for our QCD understanding of exclusive 
reactions. 
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1 Introduction 



At a sufficiently large scale exclusive QCD reactions factorize into a perturbative calculable par- 
tonic hard-scattering amplitude and universal hadron distribution amplitudes [H 121 El • The study 
of such reactions offers the possibility to directly explore non-perturbative features of hadrons at 
the amplitude level, as well as to test our understanding of the amplitude factorization. Unfortu- 
nately, exclusive reactions are still challenging to both experimentalists and theoreticians and the 
onset of the perturbative approach is controversially discussed in the literature. 

The photon-to-pion transition form factor, appearing in the two-photon amplitude of the pro- 
cess 7*(gi) 7^(92) — * vr°(P) can serve for a thorough study of the mentioned problem. In this 
process the partonic content of a meson is only probed by the electromagnetic interaction. Fur- 
thermore, since we require that the meson is produced at light-like distances, i.e., that at least one 
photon is far off-shell, this process belongs to quite a large class of two-photon processes calcula- 
ble by means of the operator product expansion (OPE) jlj. Deeply virtual Compton scattering 
(DVCS), deeply inelastic lepton-hadron scattering (DIS) and the production of various hadronic 
final states by photon-photon fusion belong to this class of processes. Such processes can be de- 
scribed by a general scattering amplitude given by the time-ordered product of two-electromagnetic 
currents sandwiched between the hadronic states. For a specific process, the generalized Bjorken 
kinematics at the light-cone can be reduced to the corresponding kinematics, while the particular 
hadron content of the process reflects itself in the non-perturbative part of the amplitude. Hence, 
the generalized hard-scattering amplitude enables us to relate predictions of different two-photon 
processes on the partonic level. 

In the leading twist- 2 approximation, the light-cone OPE approach is equivalent to the collinear 
factorization scheme P3I21E1- The transition form factor factorizes as a convolution of the hard- 
scattering amplitude T and the pion distribution amplitude 0, with respect to the momentum 
fraction x: 

F^(w,Q) = f n T(u,x,Q,fj, f )<g)^(x,fj,f), ®=/ dx, (1.1) 



where 



n2 _ (9i ~ Q2) 2 _ g| ~ ql 

~ ' ' U ~ ,2 1 „2 • 



4 q{ + qi 

In the above, the resolution scale Q 2 is large and the asymmetry parameter uj is fixed, i.e., \uj\ < 1, 
while Hf represents the factorization scale. Due to Bose symmetry the transition form factor is 
symmetric in uj. The perturbative expansion of the hard-scattering amplitude reads 



T{uj,x,Q,n f ) = ^ 



T<°> (u, x) + <^TVL x, Q) + ##T( 2 4 x, I, 9. 

2tt \ fi f J (2tt) 2 \ n f fi r 

2 



+0(ai) + {x ->l-a;} , (1.2) 
where \x r is the renormalization scale and the leading-order (LO) contribution is given by 

(L3) 

The normalization of T given above corresponds to <j>(x, ///) normalized to one and f v = 131 MeV. 
Note that a residual dependence on the renormalization scale fi r appears in the truncated pertur- 
bative expansion of the hard-scattering amplitude. The next-to-leading order (NLO) correction 
to the hard scattering amplitude has been calculated in the modified minimal subtraction (MS) 
scheme El- in the next-to-next-to-leading order (NNLO), only the contributions coming 
from the quark bubble insertions have been evaluated [HI El ED], again using the MS scheme. The 
pion distribution amplitude <ft(x,fif) is intrinsically a non-perturbative quantity and cannot be 
determined from the perturbation theory. However, its evolution is governed by the evolution 
equation 

in which the evolution kernel has a perturbative expansion as 

V(x,y^ f ) = ^V^(x,y) + ^-V^(x,y)+0(al) . (1.5) 

The evolution kernel has been estimated to NLO accuracy using the MS scheme [TTj IT2| IT3] and 
the corresponding solution of the evolution equation has been obtained E3 HO] - The latter 
can be expressed in the form 

<Kx,» f \no) = (o W/M + ^^0 (1 W/|/io) + O (a>) , (1.6) 

where the scale /j,q denotes some low scale at which the non-perturbative input has been obtained. 
The solution (jl.fi)) satisfies the initial condition <j)(x, ^q|//o) = <jft\x, ^o\^o) an d for [if — > oo takes 
the asymptotic form 0(x, \ij — > oo|/x ) = 6x(l — x). We stress that the evolution equation as 
defined by (jl.4j) and (jl.5)) corresponds to the simplified scheme fixed by the preference that the 
distribution amplitude <fi should have no residual dependence on the renormalization scale 1 . 

The photon-to-pion transition form factor has been measured at large Q 2 by the CELLO [17\ 
and CLEO JH] collaborations, where one photon is almost on-shell, while the second one has a 



1 Note that, in general, such a residual dependence appears along with the evolution kernel depending on two 
scales: 

V(x, y,M/) — ^ V® (x, y) + ^ (V« (*, y) - | V (0) (x, y) In ( + O (a») . 
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virtuality up to 9 GeV 2 . Different authors analyzed the data [H |13 IZD 122 123 12S 123 113 and it 

is often stated that the pion distribution amplitude is close to its asymptotic shape; for previous 
work see also [23 EBl EH| ■ However, in this kinematics, the shape of the distribution amplitude can 
be constrained only from the scaling violation that arises from the evolution of the distribution 
amplitude. In the small |o;|-region, perturbative QCD gives a parameter free prediction of the 
photon-to-pion transition form factor and in the intermediate region, one might extract a few 
lowest moments of the distribution amplitude and confront them with non-perturbative results, 
see Ref. j3U]- However, both high precision data as well as a precise understanding of perturbative 
and non-perturbative effects are necessary for such analysis. 

Thus, the computation of both perturbative and power suppressed contributions is an im- 
portant task. In this way, we can gain insight into the perturbative approach to exclusive pro- 
cesses. However, calculations of exclusive amplitudes beyond the leading order are quite cum- 
bersome. Besides the photon-to-meson transition form factor and similar exclusive two-photon 
processes, the perturbative next-to-leading order predictions are known only for the pion form 
factor jHIl E21 HHH HHl E3 HUH IHZj and for the amplitude of charged meson pair production in 
two-photon collision for the case of equal momenta sharing meson distribution amplitude [3H] . 
Fortunately, in the perturbative sector massless QCD is invariant under conformal transformation 
provided the coupling has a fixed point, so that the /3-function, the renormalization group coeffi- 
cient of the running coupling, vanishes. In lowest order of a s , the conformal symmetry breaking 
part, which is consequently proportional to (3/g = p a s /(4:Tr) + O(a^), can be determined by 
calculating the abelian part of the gluon self-energy proportional to the number of quarks n/. 
Additional subtleties may appear owing to the factorization procedure and they can be resolved 
by a finite renormalization of the hard-scattering and the distribution amplitude. Making use of 
conformal symmetry constraints, together with the explicit calculation of terms proportional to 
the /^-function, offers a considerable simplification of the perturbative calculation and, in our case, 
gives the possibility of going beyond the NLO approximation. 

Indeed, for the photon-to-pion transition form factor we can take advantage of this symme- 
try and its predictive power (SU EHj by means of the conformal OPE (COPE) |4*T1 I42j . in which 
the form of the Wilson-coefficients is constrained. The normalization of these coefficients can 
be recovered in the forward kinematics from the DIS results for the non-singlet coefficient func- 
tion of the polarized structure function g\ known to NNLO jlSI- This field-theoretical approach 
has been explored [TSJ EH EH E3 EE] and tested to NLO |4*7j . where the /3-function is absent 
in the Wilson-coefficients. We emphasise that the "conformal symmetry breaking" due to the 
factorization procedure in the MS scheme and the restoration of conformal symmetry by finite 
renormalization are well understood at NLO. Further consistency checks are based on comparison 
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with explicit results (e.g. hard-scattering amplitudes for two-photon processes in the light-cone 
dominated region, the flavour non-singlet kernel, quark bubble insertions in singlet kernels), and 
with constraints coming from the N = 1 super Yang-Mills theory ^H] • 

In this paper we apply the COPE combining the NNLO result for the non-singlet coefficient 
function of the polarized structure function gi with the explicit result for the tif (i.e., /?o) 
proportional NNLO contribution to the hard-scattering amplitude of the photon-to-pion transi- 
tion form factor ED] , both being evaluated in the MS scheme, to obtain a full NNLO result for 
the photon-to-pion transition form factor in the so-called conformal factorization scheme. Alter- 
natively, we propose a scheme in which the hard-scattering amplitude can already be constructed 
from a knowledge of the non-singlet coefficient function of g\ and the corresponding anomalous 
dimensions. 

The paper is organized as follows. For the convenience of the reader, in Section 2 we review the 
predictive power of conformal symmetry relevant to the photon-to-pion transition form factor. We 
then propose two treatments of terms proportional to the /3-function and discuss the remaining 
freedom in the choice of the factorization procedure. The general structure of the hard-scattering 
amplitude in the MS scheme, and the NNLO term that is proportional to /3q, are analyzed in 
Section 3. For the phenomenologically important case, in which one photon is quasi-real, we 
evaluate conformal moments for the hard-scattering amplitude and making use of the NNLO 
results for the gi function we obtain the NNLO prediction for the photon-to-pion transition form 
factor in the conformal factorization scheme. We extend this procedure to other values of photon 
virtualities and present a detailed investigation of the conformal partial wave decomposition of the 
transition form factor in different u regions. Based on these results, in Section 4 we analyze the 
size of the NLO and NNLO effects for one quasi- real photon and in the small and intermediate |u;| 
regions. Finally, a summary and conclusions are given in Section 5. The appendices are devoted 
to technical details: the Feynman-Schwinger representation of the hard-scattering amplitude, 
a consistency check at NNLO between the n/-proportional MS results for the hard-scattering 
amplitude of the photon-to-pion transition form factor and the results for the non-singlet coefficient 
function of the DIS polarized structure function gi, evaluation of the conformal moments of the 
hard scattering amplitude, the Taylor expansions in u, and the prescription for reconstructing the 
hard-scattering amplitude in the momentum fraction space from the known conformal moments. 

2 Outlining the conformal symmetry formalism 

In the physical sector, massless QCD at the tree level is invariant under conformal transforma- 
tions, i.e., under space-time transformations containing the Poincare transformations, dilatation, 
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and special conformal transformations. The latter are composed of translation, inversion, and 
translation again. Conformal symmetry implies an improvement of the energy-momentum ten- 
sor, which then becomes traceless. Owing to this symmetry, one has additional constraints for 
field-theoretical quantities, e.g. for Green functions. This subject was intensively studied in the 
60th and 70th in the four- dimensional field theory. In QCD, conformal symmetry is manifested 
in the Crewther relation and in the solution for the mixing problem of composite operators 
under renormalization [SU1E3EE]- The reduced matrix elements of the conformal operators O, 
sandwiched between the vacuum \Q) and one pion (ir(P)\ states, are pertinent to the expansion 
of the distribution amplitude 2 

<Kw) = £' X -^wf lc T^ x ~ i){<P)\o n {H)W eA i = {3+ ^fll ) 2) - (2-1) 

Here cf 2 are the Gegenbauer polynomials with the index v = 3/2 of order j and the sum runs 
over even j. In this representation, the transition form factor reads 

oo 

F^(u;,Q) = U J2' T ^Q^f) (^P)\O n ^ f )\ny ed , (2.2) 

j=0 

where Tj denotes the jth conformal moment of the hard-scattering amplitude: 

T 3 (u,Q^ f ) = £ dxT(u,x,Q,fM f ) Xil ~* ] Cf(2x-1) (2.3) 



V2 
3Q 2 



27r V v-s) ( 27r ) V v-s Vt) 



As reviewed in Section 12. 1| the operator mixing problem under renormalization beyond the 
one-loop level is solved by the restoration of conformal symmetry. In Section 12.21 this allows us 
to employ conformal symmetry to the OPE of two electromagnetic currents, and to fix the hard- 
scattering amplitude ()2.3|) in the conformal limit. Additionally, we include /3 proportional terms 
and in Section |2~31 we discuss the corresponding ambiguities of this procedure. The solution of the 
renormalization group equation to NNLO is worked out in Section 12.41 



2 



It is common in the literature that the distribution amplitude is expanded in the form 

oo 
3=0 

where Bj (j — 2,4, • • •) essentially represent the non-perturbative input. Comparing with our definition l|2.1|l . we 
have (ir(P)\Ojj(nf)\£l) icd = 6NjBj(/if), where (n (P)\(D oo {[!■ f)\Q) red = So = 1 is a renormalization group invariant 
quantity. 
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2.1 Renormalization properties of conformal operators and the con- 
formal scheme 



Let us start with the constraints for the renormalization of composite operators. In the flavour 
non-singlet and parity odd sector the twist-2 operators read jHH E21 ESj 



O fl = i>{x) (n ■ 7) 75 Cf 2 ^ ) (<n ■ <3)< (2.4) 

1 n • O x=y=0 



where d=d + d, D=D — D and n is a light-like vector that makes these operators symmetric 
and traceless. The Gegenbauer polynomials C^ 2 arise from the group-theoretical construction of 
the operators and they are of order j, where this label is related to the conformal spin j + 1, i.e., 
the eigenvalue of the Casimir operator of the so-called collinear conformal group. These operators 
have spin I + 1 and canonical dimension / + 3. In other words, we have different infinite irreducible 
representations of the conformal algebra, called towers, that are characterized by the conformal 
spin j + 1, while the members of each representation are labeled by the spin / + 1. The conformal 
operators with I = j are the lowest members of each conformal tower, and we can climb the tower 
by acting with the generator of translation. 

Employing Poincare invariance, the general form of the renormalization group equation for the 
operators introduced above reads 

^°3i = -ill3kO k i. (2.5) 
d ^ k=0 

In the conformal invariant theory operators of different conformal towers do not mix under renor- 
malization. Indeed, the anomalous dimension matrix 

7i* = + + + 0(at) , with 7iS7i „ (2.6) 

is diagonal at LO. This property is induced by conformal symmetry at tree level. The fact that 
these operators will mix beyond LO even for the vanishing /5-function in the MS scheme, has been 
considered as an unexpected breakdown of conformal symmetry. Note that the appearance of the 
anomalous dimension already requires a 'redefinition' of the conformal representation at tree level, 
i.e., the scaling dimensions of the operators change. 

The understanding of this subtlety is the key for the application of conformal symmetry in all 
orders of perturbation theory. It is well known that, on the quantum level, conformal symmetry 
is broken owing to the regularization of ultraviolet divergences, which shows up in the trace 
anomaly of the (improved) energy-momentum tensor. This trace anomaly is given as a linear 
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combination of different renormalized operators. In the dimensionally regularized theory with 
space-time dimension n = 4 — 2e, it reads 

^ (G%(x)) 2 + ■ ■ ■ , (2.7) 



where the /3-function in the regularized theory is defined as 

p e = fIj f = — <? + /3 with ^ = -±#, + 0(0$ 



(2i 



and /3q = (2/3)n/ — (11/3)C^. Besides the square of the renormalized field strength tensor G^, 
multiplied by the /3-function, the trace anomaly ()2.7|) contains equations of motion and BRST- 
exact operators. Therefore, it is sometimes believed that in the physical sector of the theory the 
breaking of conformal symmetry is in general proportional to the /^-function. However, if one deals 
with composite operators, the operator product and the trace anomaly of these operators contain 
additional ultraviolet divergences. Being multiplied by the (4 — n) contribution in /3 e , these UV 
divergences produce finite symmetry breaking terms that are not proportional to the /3-function. 
The appearance of anomalous dimensions of composite operators can be understood in this way, 
too. 

A detailed analysis shows that the non-diagonality of the anomalous dimension matrix ob- 
served in the MS scheme at NLO originates from such an effect of conformal symmetry break- 
ing. It already appears at LO in the Ward-identities of these operators with respect to the 
special conformal transformation. The calculation of this special conformal anomaly matrix 
7 C (Z) = (a s /2vr)7 c (°)(/) + 0(a 2 s ) results in 



7 c(o) (O = -K0f u, + w 



(o) 



(2.9) 



where 



and 



A 



C F { 



b jk (l) 



2(1 + k + 3)5 jk - 2{2k + 3) j - k > and even 
otherwise 



x 



A(2k + 3){j-k){j + k + 3) 
A jk -i/>(j + 2)+i/>(l) 



(k + l)(k + 2) 



+ 



2A; 



(j-k)(j + k + 3) 







j — k > and even 
otherwise 



(2.10) 
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with ip (z) = jMnr(z) and Cp = 4/3. It induces off-diagonal entries in the anomalous dimension 
matrix ([2.6)1 : 



where a(j,k) = 2(j — k)(j + k + 3). The prediction ()2.11|) coincides with the moments of the 
explicitly calculated evolution kernel at NLO [TT| IT2*1 ITS]. On the other hand, the diagonal part 
of the anomalous dimension matrix ()2.(ij) 

_ a s (0) . «s (1) , a s (2) . r,, zu /„ 10 \ 

7ii=7i = ^7] + ^7, +^7j +0(«.) (2-12) 

coincides with the anomalous dimensions of the operators that are restricted to the forward kine- 
matics and are known as the moments of the splitting kernels in DIS. 

As we can see, there is a complete understanding of conformal symmetry breaking in the MS 
scheme. Now the questions arise: 

• Can we find a scheme at which conformal symmetry holds true? 

• Can we then use the predictive power of conformal symmetry? 

The first question has a positive answer in the case that the /3-function has a fixed point. 
Instead of relying on the hypothetical fixed point, we simply freeze the coupling by hand, which 
implies (3 = 0. It is then possible to find a scheme in which the renormalized conformal operators 
(|2.4|) form an irreducible representation of the collinear conformal group, i.e. their special confor- 
mal anomaly and anomalous dimension matrices are simultaneously diagonal. The rotation from 
the MS to such a scheme, which we call the conformal subtraction (CS) scheme, is given by the 
matrix B defined by 

O cs = B^O**, B jk = 5 jk + ^B$> + 0{a]) . (2.13) 
The NLO term Bj k is entirely determined by the special conformal anomaly ()2.9j) and reads 

c(0) 

a(j, k) 

Here we introduce the notation dj k = bj k /a(j,k) and gj k = Wj k /a(j,k). In the case of a non- 
vanishing /^-function, an additional off-diagonal term appears and, thus, the complete anomalous 
dimension in the CS scheme reads [H3 



Y {0) + AA 7 (0) 



B§> = -9{j >k)^f r = B{j > k) (J 7 (°) - g] . (2.14) 



7 cs = b- 1 ^ ms B + B- 1 



d . 
H—B 

djj, 



5 jklj + 9(j>k)^-A jk . (2.15) 



jk 9 



The addenda of the anomalous dimension matrix (j2.15|) is known in the lowest order of a s [40J 

Note that the diagonal part (|2.12|) . which corresponds to the moments of the splitting kernels, 
also includes (3 proportional terms. The moments of the splitting kernels are completely known 
to two-loop accuracy |5H 157)1 EHl EZj and the first thirteen entries at three-loop level are given in 
Their scheme dependence is fixed by the fact that they are evaluated in the MS scheme. 
This implies that the diagonal part of B is given by the identity matrix. 



2.2 Conformal operator product expansion 

Let us now turn to the second question we have raised. As we have discussed in the preceding 
subsection, one advantage of the conformal scheme is that, up to off-diagonal terms proportional 
to /3, the anomalous dimensions are fixed by the DIS results, and, in the conformal limit, they are 
partly known to NNLO order 3 . Furthermore, the class of two-photon processes that are light-cone 
dominated, i.e., for Q 2 large, can be treated by means of OPE. That includes the evaluation of the 
corresponding non-forward Wilson-coefficients. Under the assumption that conformal symmetry 
holds true, these coefficients are fixed up to normalization factors that coincide with the Wilson- 
coefficients appearing in the deep inelastic scattering structure functions F\ and g\ |4T| I42j . Hence, 
in the conformal scheme, taking the conformal limit in which conformal symmetry holds true, we 
can use this predictive power of the COPE to avoid cumbersome higher-order calculations. Indeed, 
the NLO coefficient functions for the hadronic tensor in the general off-forward kinematics were 
predicted in this way jH] and they coincide, after rotation to the MS scheme, with explicitly 
calculated ones jHOJ EU E2] • 

For the photon-to-pion transition form factor, the leading twist-2 result of the OPE is given 
by Eq. (|2.2|) . where the conformal moments Tj of the hard-scattering amplitude correspond to the 
Wilson-coefficients Cj, which are conventionally normalized as 

Tj(u, Q, fi) = —Cj{uj\a s {n), Q/fi) . (2.17) 

As we have mentioned, in the formal conformal limit the Wilson-coefficients are constrained in 
the CS scheme by the predictive power of the COPE: 

a s (fi) a* — fixed implies the reduction Cj(u\a s (n), Q/fi) =>- C j{oj\a*,Q / ji) , 

3 In principle, we then also know the Efremov-Radyushkin-Brodsky-Lepage (ER-BL) evolution kernels, which 
can be obtained from the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) kernels through an integral trans- 
formation |59j . 
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where Cj(u\a*,Q//j,) is given by 

. 1 . . . 

2r] uj 



^-^"•Ho'J (i+^y«^ B(j + 1 ' J+2)2Fl l 20 + 2 + H) 



1+770; 



(2.18) 



In the limit 77 — *■ 1, one obtains the result for the production of a (pseudo) scalar meson by two 
virtual photons, while for 77 = the forward case is reproduced 4 . Note that Cj{uj) is an even 
(odd) function of uj for even (odd) j, which is guaranteed by the linear transformation properties 
of the hypergeometric functions 2^1- The normalization Cj(a s ) coincides with the flavour non- 
singlet Wilson-coefficients of the polarized structure function gi taken at fi = Q. It is given as a 
perturbative expansion 



and is known to NNLO in the MS scheme jinj. Strictly speaking, this coincidence appears just at 
the hypothetical fixed point a s = a*. However, since we know the forward anomalous dimensions 
and the Wilson-coefficients perturbatively in the MS scheme, we can easily restore (3 proportional 
terms in these quantities beyond the conformal limit. 

Conformal symmetry breaking terms, proportional to the /3-function alter the COPE result 
()2.18|) in the full theory. Obviously, in the irreducible conformal representation used, the f3 term 
cannot be fixed from the requirement of conformal invariance. Thus, the definition of the conformal 
scheme 

C cs (a a (n), Q/fi) = C m (a s (fi),Q/ft) B (a.Qi)) with C cs (a* s , Q/» = C (a* s , Q/» , (2.20) 

is ambiguous and, consequently, the I?-matrix is uniquely defined only up to (3 proportional terms 

that are off-diagonal. At the two-loop level, C contains /3q proportional terms appearing in both 

(2) (2) 

Cj and 7] . Let us first set such terms to zero and, with the help of this, single out all (3q 
proportional terms in C cs : 

C cs = C m B\ p=d0 + ^8C (2.21) 



9 

C\ P=0 +-SC, with 6C=£(C m B-C m B\p-_ (h 



9 ' P 

and a s remains running. At LO a f3$ term does not appear in C MS , thus, the perturbative 
expansion 

SCj = 5Cf ] + ^5Cf ] + 0{a 2 s ) (2.22) 



4 The Wilson-coefficients appearing in the deep inelastic scattering structure functions F-y and g\ are derived 
in the usual D1S operator basis, which differs in the normalization from the definition for the basis of conformal 
operators. Hence, these Wilson-coefficients differ slightly from the r\ = limit of l|2.18|l . 
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holds true. Note that since B$ = 5 jk (l2~T3jl the LO coefficient C MS(0) = C (0) is independent of 
the scheme. Similarly to (j2.21j) . we can write the matrix B in the general form 

B = B\b= q + -6B with 5B = 5B {0) + ^5B W + 0(a 2 s ) . (2.23) 
g Lix 

If we define B so that it contains no pVterm in order a s , i.e., if we take 5B^ = as in our 
definition of the CS scheme (j2.13H2.lH) . the coefficients 5C up to NNLO read 

5C^ = 0, SO® = |- (C m V> - C m (% 0=0 ) + C<®5BM. (2.24) 

Po v J 

Since we have required that the diagonal entries of B should be one, the normalization coefficients 

Cj coincide in the forward limit with the Wilson-coefficients of g\ calculated in the MS-scheme. 

2.3 Ambiguities in the definition of the conformal scheme 

As we have discussed, the ambiguity left in the definition of the conformal scheme in the full theory 
resides in the (3 proportional off-diagonal terms, i.e., in the choice of 6B in (J2.23J) . Adopting the 
definitions ()2.13|) and (|2.23|h we set SB^ = and in the following discuss different choices of 
5B^\ restricting ourselves to NNLO. 

2.3.1 Defining CS and CS schemes 

The naive choice is to set 

5B (1) = 0, (2.25) 

which means that the (3q proportional term in the NNLO Wilson-coefficients, i.e., 5C^\ is entirely 
evaluated in the MS scheme. Since the conformal symmetry breaking part appearing in the COPE 
has to be proportional to the /3-function, we can calculate SC^ by evaluating the contributions 
proportional to the rif piece of po, i.e., from the one-loop Feynman graphs with an additional 
quark-bubble insertion. 

In this naive scheme (12.25)1 . which we denote by CS, both the conformal operators and the 
Wilson-coefficients will mix under renormalization to NLO accuracy owing to the running of the 
coupling. Let us consider this in more detail. Since the transition form factor is invariant under 
renormalization, i.e., 

fi— F^(u, Q) = 0, (2.26) 
an 

these effects will compensate one another. Thus, the renormalization group equation for the 
Wilson-coefficients reads 



d B 
an ag 



OS 



(V) ogQQ (i)' 
2vr lj (2vr) 2 Jj 



C cs + P g/^0 cf s A g) + o(a 3 )j (2 . 27) 

9 i=j+2 111 
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where the addenda is defined in Eq. (J2.16|) . Since Cf s is conformally covariant to the NLO 



i j 



approximation, i.e., it contains no partial waves (J2.18j) with a conformal spin larger than j + 1, the 
off-diagonal entries on the r.h.s. arise from the explicit fx dependence of C^ 2 \ which has been 
taken from the MS scheme. 

In the alternative conformal scheme, denoted in the following by CS, this intermediate mixing 
is avoided by the complete diagonalization of the renormalization group equation. This can be 
achieved by including an explicit fx dependence in the 5-matrix, i.e., by taking 



,*2 ' 



ln vw ) A % 9{J >k)+ 



(2.28) 



In the order we are considering, the matrix B now reads 



B 



Sjk + 



2tt 



(2vr) 2 



R (2) , 00 

B jk + - 



*2 s 



In 1^1 Af>(j > k) + 



+ 0(a 3 s ). (2.29) 



with Af fc } 



A CS(0) 



0, i.e., LO addenda 



This choice introduces 7 CS given by Eqs. 
to the anomalous dimension matrix vanishes in this scheme. However, for dimensional reasons 
the choice ()2.29j) additionally introduces a new residual scale dependence fi*. The meaning of 
this procedure is obvious. We do not resum the remaining off-diagonal ln fx terms through the 
renormalization group equation, rather we include them in the Wilson-coefficients, where they 
will be annulled. This is indeed on the same footing with that we have already discussed in 
the [3 — case, where the off-diagonal entries present in the MS scheme have been removed 
by a finite, however, fx independent renormalization. Consequently, both the operators and the 
Wilson-coefficients 



CS = 



c 



cs 



C MS B, 



(2.30) 



where B is defined by Eq. (12.29)) . now satisfy the desired renormalization group equations in the 
a 2 s approximation: 



a s (fx) 



(2.31) 



dfx dg 



C 7 CS (a s (fx),Q/fx) 



2tt 



-7 ' + 



(27T) 



Cf s (a s (fx),Q/ fx) (2.32) 



Note that the forward anomalous dimensions jj remain explicitly ft independent. However, both 
the off-diagonal piece of the anomalous dimensions and the Wilson-coefficients now possess a 
residual fx* dependence at the orders and a 2 s1 respectively. 

To restore the fx dependence of the Wilson-coefficient, we perturbatively solve its renormaliza- 
tion group equation (|2.32J) . Up to an integration constant 5C, its solution to two-loop accuracy 
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can be expressed by the Wilson-coefficient ([2.18)1 . appearing in the COPE, 

-5C' 



a 2 fa 



Cy b («.(//), Q/fi) = Cj (a a (jj), Q/fi) + 2 
which now depend on the running coupling: 



^2.33) 



Q _d_\ I // 



2 \ 2(2u>YB(j + + 2)^ /j + 1 + j + 2 + | 7i 



:i + ^y +i+ - 



2cj 



1 + 



(2.34) 



The running of a s (/i), appearing in the lowest approximation of both 7 j and Cj, is compensated 
by introducing a logarithmic change of the /x dependence in the normalization factors: 



( ( \ Q 9 

Cj \ a s /i , — ,— 



c (o) + Q«(M) c (i) + ««(m) c (2) + a 8 (/x) ^o ln ^ 




(27T 



,2 J 




5 



<9 7 



+ 0(al 



(2.35) 



Equation (|2.34j) should be understood in the sense of a consequent expansion with respect to a s 
up to the order a 2 s . 

So far we have found a rather natural way to include the effects of the running coupling in 
the structure of the COPE result, with the advantage that the conformal operators do not mix 
under renormalization in NLO. It remains to fix the integration constant SC, which vanishes 
in the kinematical forward limit. We can identify it with the non-covariant part calculated in 
the MS scheme, in an analogous way to our discussion in the case of the CS scheme. On the 
other hand, it is rather appealing that the Wilson-coefficients Cj (a s (fi),Q 2 / fi 2 ) contain only 
conformally covariant terms to NNLO. For the scheme we call CS, we adopt this prescription, 
i.e., we put 5C = 0. In the NNLO approximation, we then have a partial wave decomposition 
of the transition form factor with respect to the 'good' quantum number - conformal spin. This 
in principle allows us to extract the conformal moments of the distribution amplitude with a 
'well-defined' conformal spin for experimental data. 



2.3.2 Calculational prescriptions 

Let us comment on the renormalization scale dependence and give the calculational prescription 
for the schemes we have proposed. 

First, we have introduced a naive recipe (CS scheme) which combines the COPE result with 
that one explicitly calculated in the MS scheme. Now we extend our analysis by distinguishing 
between the renormalization scale /i r (the argument of a s in the Wilson-coefficients), and the 
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factorization scale fj,f. We require that the matrix elements of conformal operators should depend 
only on the factorization scale fif. Thus, the scheme transformation now reads 



C cs (a s (fi r ),Q/fi f ,Q/fi r ) = C Mti (a s (fi r ),Q/fi f ,Q/fi r ) B(a.(jif)). 



Employing the scale-changing relation 

a a (jj, f ) = a s (fi r ) 
we expand the rotation matrix 



a 8 (fi r ) Pq 

1 H in 

2vr 2 



A, 



B(a s ( H )) = B(aM) + ^^\n 



0(ai 



B (1) +0( 



Hence, in this scheme, the Wilson-coefficients read to NNLO accuracy as 
C cs = C(a s (ti r ),Q/nf) 



A) atl(jir 
2 (2vr)< 



_ C rMS(2) | 



(2.36) 



(2.37) 



(2.38) 



a s (fi r ),Q/fi f ,Q/fi r ) +ln ( ^ 1 C^B^ 



+ 0(a 3 s ), 



where C^f S ^ denotes the {—(3q/2) proportional contribution evaluated in the MS scheme, while 
the structure of C is fixed by (j2.18H2.l9j) . 

Alternatively in the CS scheme, we have employed renormalization group invariance to incor- 
porate the running of the coupling into the generic structure of the COPE result and have used a 
finite renormalization to preserve the structure of the COPE to NNLO accuracy: 



Q/v) = Cj (a,(//), Q/fi) 



(2.39) 



where Cj is defined by Eqs. (l2.34H2.3Kj) . The form of the Wilson-coefficients in which the distinc- 
tion between the scales \i r and \if is introduced can be obtained analogously to the previously 
discussed case of the CS scheme, and will be presented in Section I3~3l 

In Section H] we employ both of these schemes to estimate the size of NNLO effects at a 
given input scale. The missing ingredient for a consistent NNLO analysis including the evolution 
of distribution amplitude, is the anomalous dimension matrix at three-loop level. Whereas the 
first few diagonal entries have been calculated, the off-diagonal part is unknown. It could be 
read off from the n/-proportional part of the special conformal anomaly matrix in the two-loop 
approximation. Also, generally, the trace anomaly will affect the COPE at NNLO accuracy, see 
Eq. (j2.7j) . We rather make use of the freedom to rotate the conformal symmetry breaking piece 
from the perturbative sector to the non-perturbative one, as has been done in the CS scheme, or 
the reverse. However, we expect that the mixing effect in the CS scheme will be negligibly small 
and its detailed investigation is beyond the scope of this paper. 
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2.4 Evolutional behaviour of conformal operators 

We end this section with a short review of the evolutional behaviour of the conformal operators 
from which that of the distribution amplitude 4>(x, fi) can easily be established. For the convenience 
of the reader, we repeat here the basic steps for solving the renormalization group equation (|1.4jl 
and present the results in a form convenient for the phenomenological analysis fSJ HH] . 

The renormalization group equation ()2.5|) is an inhomogeneous first-order partial differential 
equation and after sandwiching the conformal operators between the hadronic states of interest 
we obtain the evolution equation for the reduced matrix elements: 

A^-^<^(^)l<^^-(A^)l^> reci = — /^C/^X^C^)!^^^/^)!^)^ — ^- fc (/x)<^(^)|cr^ A;fc ( Ai )|^>— ^. (2.40) 



The solution can be achieved by the ansatz 

(7r(P)|O ii ( A *)|n>" d = E'M/^expj- ^ lh (A (rr(P)\O kk ^ )\Qy ed , (2.41) 

i . r> J On U 



k=0 

with the initial condition 



Bjk(n = A*o, /x ) = 8jk ■ (2-42) 

The recursive solution of this set of differential equations, starting with the homogeneous one for 
j — 0, has been written for an arbitrary scheme in a compact form (see Ref. |15j): 

B» = s = S jk + £ 7 jl D + C ( 7 ND £7 ND ) Jk + " " " , (2-43) 

where 7 ND represents the triangular off-diagonal matrix 

7 nd = ( 7# for j > k p 44) 

\ otherwise 

and the operator C is an integral operator whose action is defined by 

£>ifk = -T d 4 7 £V) exp {- f % bW) - 7*(AOl} ■ 



(2.45) 



In the MS or CS scheme, the anomalous dimensions (J2.6)) depend only implicitly on the scale \i 
via the running coupling a s (fi). For the /3-function, we employ the expansion [03] 

Po+ (A /, ^ A; + 0{a s ), (2.46) 



5 4tt (4vr) 2 ^ i (4tt) 3 

2 „ 38 325 2 5033 2857 

A) = 3^/ - 11, A = y% - 102, & = — + -jg-rif - — 
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Since the off-diagonal entries of the anomalous dimensions give only subleading log's, which will 
not be resummed, we expand the B matrix in powers of a s : 



Performing the integrations in Eq. (|2.43jl . we obtain the desired results for 



^2.47) 



B jk = -#ifcO,Mo|l 



ND(1) 

Po 



[2.48) 



(2) 



(1) (1) a (°) (°) 
Tj - 7fc Pi Tj ~ 7fc 



Po 2Po Po 



ND(1) 



Po 



+R jk ((jL, (jlq\2) 



„ ND(1) ND(2) 
Pi Ijk Jjk 



2Po Po Po 



j—2 ND(1) , . , -. . , I , . 

+ yy Tjm it mfc (/X,/i |l) - R jm (fi,fl Q \2) 7, 



ND(1) 
mk 



m=k+2 



Po 



1 + 



7 7 (0) -27i 0) +7[ 0) 
ft 



where 



Po 



a i (o) (o) 



« s (/io) 

a s (/i) 



00 



(2.49) 



The leading log's, associated with the diagonal entries, are resummed, while the subleading ones 
are expanded with respect to a s : 





7 (°) 

lis 




00 







where 



Ak\^^o) 



_ a s {fi ) 



2tt 



a (0) (1) 
Pi Ik Ik 



1 + ^^Wo) + ^Af^+OW 



^2.50) 



1 r 



2Po Po Po 



(2.51) 



PI ~ PiPo if Pi lt } , 7^' 



(!) .,( 2 ) 



+ 



8/3 A) 4/5 /3 2/3 
To the considered order, the evolution of the matrix elements then reads 



(7r(P)|O ii (A*)|n> wd = E' 

fc=0 



«s(/f) 

« s (a*o) 



(0) 



a s (jj) 



2vr 



(2vr) 2 



^2.52) 



x 



red 



The off-diagonal entries are known only at NLO and are given in the MS scheme in Eq. (j2.11j) . 
In the CS scheme, they are proportional to Po, as given in Eq. (|2.15|) . whereas in the CS scheme 
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they are equal to zero by definition. Therefore, the mixing of operators is an a 2 s suppressed effect: 



B 



(i) 

jk 



(2.53) 



(o) 



jk T? 



(0) 



/3o + 7f 




— R jk (fi, ^ \2) 



1 



- J R ife (/i,/i |2)A 



CS(1) 
jk 



where A^ is defined in Eq. (I2.16)) . Here we have taken into account the explicit /^-dependence 



in the anomalous dimensions, induced by the transformation (J2.29)) . The addenda A^ S( 1 



jk 



27jfc / Poi which is presently unknown, has to be evaluated at /i = /i*. Note that the a s -power 
counting remains correct as long as the scales //*, fi, and /io are of the order of Q 3> Aqcd- The 
auxiliary scale can now be set, e.g., to fi* = \i. 

Let us remark that the evolution of the distribution amplitude can be formally obtained by 
resummation of the conformal partial waves given in Eq. ()2.1|) . Taking into account the evolution 
of the reduced matrix element in Eq. (j2.52j) . one finds the eigenfunctions of the evolution equation, 
expanded with respect to the Gegenbauer polynomials [To] . 



3 Hard-scattering amplitude to NNLO accuracy 

In the preceding section we have outlined the structure of the conformal predictions in the con- 
formal momentum space and in this one we analyze the structure of the NNLO results in the 
momentum fraction representation. In Section 13.11 we derive a convolution representation of the 
NNLO term proportional to (3q and also give the general structure of the hard-scattering am- 
plitude in the MS scheme up to the NNLO order. Furthermore, in Section 13.21 we consider the 
phenomenologically important case of the asymmetry parameter |u;| equal to 1 and in Section 
13.31 we then present the NNLO result for the conformal moments in the CS and CS schemes. In 
Section EH we analogously present the results at small and the intermediate values of \u\. 

3.1 /^-proportional NNLO terms and the general structure of the 
NNLO results in the MS scheme 

First, we consider the term proportional to rif, i.e., /3q, appearing in the NNLO calculation of the 
two-photon hard-scattering amplitude (jl.2|) . For the case of general Bjorken kinematics, the result 
has been given in Ref. [Oj and is easily restricted to the kinematics of a particular process, i.e., in 
our case, to the kinematics of the photon-to-pion transition form factor (see Appendix [B] for the 
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definitions of generalized Bjorken kinematics). In the special case |a>| = 1 (one photon on-shell), 
these results coincide with the results from Ref. jlOj . 

The authors of Ref. [2] have presented the regularized results in terms of hypergeometric 
functions 2^1- It is instructive to rewrite it as a convolution of the amplitude 

1 



T^(u,x\e) 



(1 - co(2x - 



(3.1) 



with the kernels 



v a (x,y\e) = 9(y-x) 



x 



l+e 



+ 



,yj [y-^y) 

' x\ 1+e 1 (x^x 



v (x,y\e) = 8(y - x) 

\yj y~ x [y^y) 

1 ( x\ ^"^ e ~^ cr / \ ( X — ► X 

g{x,y\e,a) = 6{y-x)-[l--} B ~, 1 + a, -e - a }+ { 

y \ y \y J [y^y 



(3.2) 
(3.3) 
(3.4) 



where B(x, a, b) = Jq dy y a ~ l (l — y) h ~ l is the incomplete Beta-function. Here a and e are the 
dimensional regularization parameters in = 4 — 2er[e]) associated with the quark bubble insertion 
and the overall loop, respectively. The kernels [v a (x, y\e)] + and v b (x, y\e) ^ are diagonal with 
respect to the Gegenbauer-polynomials C^ 2+e (2x — 1) and they are regularized with the usual 
[•••]+ prescription 



[v(x,y)],=v(x,y)-6(x-y) / dzv(z,y) 

Jo 



Their eigenvalues are 



l + e 



1 



v b i = 2ip(2 + e)-2ip(2 + e + j). 



(3.5) 



(3.6) 



J (l+j + e)(2+j + e) 2 + e' 

The g kernel is not diagonal with respect to the Gegenbauer-polynomials and is responsible 
for the apparent breaking of conformal symmetry in the MS-scheme. Its expansion 

g(x, y\e, a) = g(x, y) + g'(x, y) e + g(x, y) a + O (e 2 , cr 2 , ea) 



reads 



g{x,y) 
g'(x,y) 
g(x,y) 



0(y 



y-x 



2 
Li 2 



+ 



x 

y 



! , 
y J 



ln z 1 - 



x 



y-x 

Q(y ~ x) 

y-x 



+ 



x 

y 



1 x 
y, 



Li 2 (i: 



■> X 

* V 



x 

y 



(3.7) 



X 



+ g'(x,y). 
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There is a similar expansion of the v 1 kernels (i = a,b) 



v'(x, y\e) = v\x, y) + v*(x, y) e + -v\x, y) e 2 + O (e 3 



where 



v*(x,y) = 9{y-x)f i (x,y) + 



y -»■ V 



v\x, y) = 8(y - x)f(x, y) In ( - j + J X ~" X \ , 

\y [y^y\ 



(3. 



rf \ I /■>' i />"■ 



^(ar, y) = 0(3/ - x)f (x, y) In 2 ( - 1 + 

\yj [y^y) 

and the functions /' can be read off from Eqs. (|3.2|) . (|3.3|) . 

The LO kernel of Eq. ()1.5|l is expressed in terms of the v l kernels introduced above: 

V (0) (x,y) = C F [v(x,y)] + , v(x,y) = v a (x, y) + v b (x, y). 

For the NLO kernel we use the colour decomposition 



V (1 \x,y) = C, 



C F v F (x, y) - -j-v p {x, y) - (c F - -£\ v G (x, y) 



(3.9) 



(3.10) 



+ 



(3.11) 



In the following we particularly need the /^-proportional kernel 

5 

vp(x, y) = v(x, y) + -v(x, y) + v a (x, y). 

The unrenormalized NLO and NNLO corrections to the hard-scattering amplitude are of the 
form 

fW{u,x) = C F fP(u,x), 



(3.12) 



f( 2 \u,x) = d 



C F fP^x) - %Tf\u,x) - (C F - ^) t f{u,x) 

Employing the integral representation of the hypergeometric functions, one can express the regu- 
larized results of Ref. [5| 

if\u,x) = M(u,x\e,0) , (3.13) 
?f\u,x) = -3 ^^^- M{u;,x\e,e) (3.14) 

in terms of convolutions 



M(cu, x\e, a) 

r(e + a)T(2 - e)T(l - e - a) /4 W +<T rl 



x < - 



r(3 - 2e - a)V(l + a) 
(l-e)(l-4e-3a) 



Q 2 



dyT(°\LU,y\e + 



(J 



(1 - e)[2 - e + 2e(e + a)] - eo(e + a 
1 -e 



x (l-e)[2-2e(l-e)+a(l-a)] a . . . 
5(x -y) + ± ^ K — '- i ^ v a (y, x\a) 



l + a 



v b (y,x\a) , -{e + a)[g(y,x\e,a)) + ) .(3.15) 
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The results given above contain UV and collinear singularities, which are removed by renor- 
malization (introduces the scale fi r ) and factorization (at the scale ///) of collinear singularities. 
The renormalization procedure in the MS scheme (for details see Ref. ^Oj) induces the following 
general structure of NLO and NNLO corrections present in expansion (J1.2j) : 



TM(u,x,Q/n f ) = C F T F l \u,x)+\n(%\ \T^®V^\ (u,x) 



,Q//i f ,Q//i r ) = C F 



r T (2) #W2) 



where 



T F \u, x) 



T f 2 \uj,x) 



if fax) 



Tg\u,x) 



2 -P 
x J T (o)^y(D + 



C F - ^] TP 



{uj,x) +m\— T 



Cf tP + - In ( ^_ ) T(°) ® 



(3.16) 
(3.17) 



/So 



'Q 2 



T<°> (a;, y) <g> {Ti 1} (y, x) + LN(w, y) [ufo, x)] + } , 

r«V,j/) <g> |ri 2) (y,x) + LN(w,y) ([M^)]+ + T« ® [v]+(y,x] 

+^LN 2 (^2/)M+® [v]+{y,x)\, 



(3.18) 
(3.19) 



T(°\u,y)®\TP(y,x) + LN{u,y) 



-r^){y,x) 



(3.20) 



-^LNV,y)Ky,x)] + j, 



T (0) (u, y) ® Tg> (y, x) + LN(u, y) [v G (y, x)] 



+ f> 



(3.21; 



while LN(u>, x) = ln[l + u — 2xuf\ and T^(lj, x) is given by Eq. ()1.3)) . For a detailed discussion of 
the appearance of the LN terms, see Appendix 1X1 

The explicit calculation provides us with the kernels 



rP{x, v ) 



—v° + g 
2 y . 



29 
— 'i 
12 



— w a + v a 



(x,y) - -5(x-y), 
f- * 

209 7.1 .. 
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V V V H (7 + (7 

36 3 4 6 



(x,y) - 3S(x - y). 



(3.22) 
(3.23) 



As we see, these T-kernels are built up of v and g kernels appearing in the evolution kernel, too. 
While the g and v kernels appear at NLO, the g and v ones show up for the first time at NNLO. For 
the missing two entries, namely, T F 2 ^ and Tq\ we expect a similar structure, but with additional 
and unknown building blocks that are related to the ^-kernel. 
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Making use of the fact that both the photon-to-pion transition form factor and the forward 
Compton scattering belong to the class of light-cone dominated two-photon processes, which can 
be described by a general scattering amplitude, we have performed a consistency check between 
the previously presented results for the hard-scattering amplitude of the photon-to-pion transition 
form factor (known up to /^-proportional NNLO terms) and the corresponding results for the 
non-singlet coefficient function of the DIS polarized structure function gi [13] • The procedure is 
presented in detail in Appendix [BJ 

A few comments regarding the LN(u;,x) terms are in order. In NLO, we observe that the 
LN(a;,a;) term match the ln-term indicated in Eq. (|3.16|h i.e., we can absorb it in Eq. (j3.18|) by 
an appropriate choice of the scale: 

-> J% = /i 2 (1 + ^ - 2XCU)- 1 with i = {f, r}. (3.24) 

The explicit NNLO result for the f3 proportional terms satisfies the same rule for the scale re- 
definition, which indicates a general property of the hard-scattering amplitude evaluated in the 
MS scheme. This is shown in Appendix El The terms proportional to LN(u,x) are vanishing in 
the limit \u\ — > and for \u\ — > 1 provide a logarithmic enhancement in the end-point region. 
However, a resummation of such terms through an appropriate scale setting is misleading, since 
other logarithmically enhanced terms also appear. For instance, at NLO we have 

\n(l-y) r . ln 2 (l-a;) + 21n(l-x) „ „ 

-J y -l ® [v(y, x)} + = —± ^ 1 + O In 1 - x)) , (3.25) 

1 — y 1 — x 

while the contribution of Tp\x,y) is [see Eq. (13.22}) ] 

.(i) ^ _ ln 2 (l -x) +31n(l -x) + 9 



® (y, x) = i '— ^ '- + O (HI - x)) , (3.26) 

l-y 2(1 -x) 

so that only a partial cancellation of the ln 2 (l —x)/(l—x) term appears. 



3.2 Limit \u\ — > 1 and corresponding conformal moments in the MS 
scheme 

Of special interest is the limit \u\ — > 1, since different (pseudo-scalar) meson-to-photon transition 
form factors are measured in this kinematical region. We can trivially perform this limit in Eqs. 
(jl.3j) . (j3.18j) . and (J3.2()|) and after convolution we present the result in the form of Ref. [TU] : 

r<°>(x) = — 1 . , (3.27) 



2(1 -x) 
1 

2(1 -x] 



ln(l - x) + - In (1 - x) 

2 2x V J 2 K ' 
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(3.28) 



2(1 -x) 



457 


/47 




_ ~ "48" ~ 


V36 ~ 


Ax J 



x) + 



13 _ 1 
12 ~ 4x 



ln 2 (l - a;) 



(3.29) 



-J ln 3 (l - x) - ^Li 2 (x) + ~Li 3 (:r) - $ 12 (x) 
o 3 2 



where the polylogarithms are defined by 
Li n (x) = 



dy 1 ^ , with Lii(x) 



2/ 

o / \ l r . ln 2 (i-y) 

Si 2 (x) = - / 

y 



ln(l — x) , 



2 Jo 

- ln 2 (l - x) ln(x) + ln(l - x)Li 2 (l - x) - Li 3 (l - x) + C(3) . 
2 



(3.30) 



Here we introduce for the special case \u\ = 1 a more convenient notation in which the ln(2) terms 
arising from the LN(a> = 1, x) functions present in Eqs. (j3.18M3.20j) are absorbed in ln(2Q 2 //x^ r - ) ) = 
ln(— qf/[if, r A We do not list here the terms proportional to these log's but refer to Ref. (TOj for 
their explicit expressions. 

The conformal moments of the individual terms we encounter in Eqs. fl3.27M3.2Uj) are given in 
Appendix O Here we list the conformal moments of T^(x), Tp\x), and Tg (x) calculated in 
the MS scheme: 



,(o) 



n(l) 



T (2) 



2j + 3 



(j + l)(j + 2) 

2J + 3 
(j + l)(j + 2) 

2J + 3 
(l+m + j) 

1 



1 



9 3-4&C7 + 1) 

2 + 2(j + l)(j + 2) + (j + l)2( J + 2)2 

1 



+ 2S!(j + l) 



(3.31) 
(3.32) 



A JL _ I C(3 ) _ Z C (2) + - 
48 2 U J 3 U J (j + i)(j + 2) V 36 



^ - ^C(2) 



+ 1) - 2sf(j + 1) - i-iys-M + 1) 



l-SfrCy + l) 
2tf + l)=»(j + 2)* 



(3.33) 
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10 

y 



' TT'^'j I 1 > ' ™< s f(j + 1) + |-S? (J +1) + i^./ I ; 



2 

3' 



2(j + l)3( J + 2)3 9 

After inspection of Eqs. (j3.1(ij) and (j3.17j) one notes that the conformal moments of the terms 
proportional to ln(2Q 2 //x 2 / can be conveniently expressed using the conformal moments of the 
kernels t> and Vp. For the definition of conformal moments of the kernels, we refer to Appendix 
O Eq. flC.7j) . The conformal moments of the diagonal kernel [v(x, y))+, given in Eq. (j3.9j) . are 
denoted by Vj = Vjf. 

3 1 



■2S& + 1) 



2 (l +J )(2 + j)' 
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(3.34) 



The conformal moments of <g> [v] + , determined using ()C.10|) . are given by T® Vj. On the other 
hand, the kernel [vp] + in Eq. (|3.11j) is non-diagonal, and as shown in (jC.8|) . both the diagonal 
vpj = vpjj as well as the non-diagonal vp^j (k > j and k — j even) conformal moments contribute 
to the conformal moments of ® [vp] + , i.e. one obtains Y,k>jT® vp^j- After performing the 
convolution and making use of the results for the conformal moments summarized in Appendix 
IU1 one can express the conformal moment of (g) [v@]+ in a form T® v^j, where we introduce 
the "effective" conformal moment of the Vp kernel amounting to 

^ = ^. +C(2) _9__J__. (3.35) 

Finally, we summarize our MS results. The LO contribution is given in Eq. ()3.31j) . the NLO 
contribution takes the form 



while the (— (3q/2) proportional NNLO term is given by 



-13,3 ^ 111 „2 



$ J + 1)0' + 2) 



^i + 2 ln (jr| ^ 



ln(f)^- (3.37) 



Note that owing to the fact that vp is non-diagonal, even the lowest partial wave, i.e., j = 0, of 
the NNLO correction depends on the factorization scale as well as the renormalization one. 



3.3 NNLO result in the CS and CS schemes in the limit \cu\ — ► 1 

Let us now turn to the conformal schemes CS and CS. We make a distinction between the 
renormalization and factorization scale. Consequently, the argument of the coupling in the Wilson- 
coefficients depends on /i r and, as discussed in Section 12. 3\ we require that the matrix elements 
of conformal operators depend only on the scale /i/. We use the COPE, where for \ou\ — 1 the 
Wilson-coefficients ()2.18j) simplify to 

C.-c.^)^— J r(i + 2 + 7 . /2)r(i + 3 + 7 . /2)r(2i + 3) - (3.38) 
The anomalous dimensions are given by 

7j = ^) 7f (3.39) 



-Up Vj + L-F 



C F v F ,j - Y r '-> - (°f ~ \°a ) 'V;., 



+ 0(al 



{ 2tt ' J ' (2tt) 

where Vj, vfj, vqj, and vpj are the diagonal conformal moments of the evolution kernels [v(x, y)] + , 
[vp(x, y)]+, [vg(%, y)]+ and [vp(x, y)]+, respectively, and they coincide with the moments of the 
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DGLAP kernels. The LO moments are given by Eq. (|3.34j) . while other entries can be found in 
Refs. jHU EH1 EI] • Analogously, we decompose the normalization factor 



a-.s 



1 + =?-C F cf + 



2^ 



0o 



(2) 



ft C, 



.(2) 



+ O (c£) . (3.40) 



Its NLO contribution reads 



C) ' = bill + 7 + -Oi 7 +2 ; — 

3 1V J> 2 u ; 2 2(7 + l)(j + 2) 



(3.41) 



while the NNLO contributions can be determined from the Mellin moments of the coefficient func- 
tions calculated from the a 2 s corrections to the polarized structure function gi [13]. A consistency 
check of the /3q proportional part of these results is given in Appendix iBl 

As discussed in Section ESI the Wilson-coefficients in the CS scheme are obtained in a straight- 
forward manner by means of Eq. (JOED, where C^ S(2) (a s (/i r ) , Q/fi f , Q/fj, r ) is given by the expres- 
sion ()3.37|) . Taking into account the proper normalization, i.e. identifying TW with by 
(|2.17j) . the expansion of Cj, defined by Eq. (|3.38jl . leads to the complete NNLO result for the 
hard-scattering amplitude: 



3Q 2 



T (o ) + a £ pi cs(i )(g/ ) + 



(2vr) 2 



Tf^ 2 \Q/ H ,Q/^) + 0{a: 



where 



T- 

3 



■CS(1) 



C F + 



2j + 3 



+ 1)0 + 2) 



In 



'2Q* 



, (3.42) 



(3.43) 



T 



CS(2) 



C F { CfT, 



Fl F,j 



CS(2) Aw2) 



2 ^ 



c F -Uj A ) rg (2) + h, 



'2Q 5 



2j + 3 



H) J (j + 1)(j + 2) 
Cf^fj + cfv 3 + tjJ (^(j + 1) + S x {j + 2) - 25i(2j + 3)) 

2Q 2 \ \ A 



+-Un 
2 



-"0 cs,s 



i ~2 C a) 'V;., 



CS(1) (Q//^/)-f ln 2 ^ g " 



C 
2j + 3 



and 



rpCS(l) 
1 F,j 



T, 



CS(2) 



aij; 0+1)0+2) j /' 

{cf + ^ [Si(j + 2) + 5x0' + 1) - 2S'i(2j + 3)]} 
(j + + 2) + (cf v s + ^) + 1) + S 1 (j + 2) - 2^(2j + 3)] 



2j + 3 

0' + 1)0' + 2) 
2J + 3 



(3.44) 



(3.45) 



CS(2) 



G,j 



+j{[Si(j + 1) + ftO" + 2) - 2Sx(2j + 3)] 2 
+S 2 (j + 1) + S 2 (j + 2) - 4S 2 (2j + 3) + 2C(2))} 
11 + " -{cgj, + woj [Si(j + 1) + + 2) - 25' 1 (2j + 3)] } 



+ 1)0' + 2) 



(3.46) 



(3.47) 
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while T^ 2 j corresponds to the MS result given by Eq. (J3.33|) . Note that, in accordance with 
the renormalization group invariance, the off-diagonal part of the anomalous dimension matrix 
proportional to /3q has been changed and, consequently, also the conformal moments 



= v,,cs,s_ ? .s (j + l)(j + 2) / CS(1) 



T^f > - T«) . (3.48) 



Finally, we present the result for the CS scheme, in which the conformal covariance of the partial 
wave decomposition is preserved. The modification concerns only the terms proportional to (3q in 
which the off-diagonal entries in Eq. ()3.44|) are removed by making the following replacements: 

CS.E 

tfS =► T TP = ( J+ 2 i)|/ +2 ) { c ^ + v n [SiiJ + + SiiJ +2) ~ 2Si{2j + 3)] } ( 3 - 49 ) 

+ ln(2) +T®±\n(2)v j . 

The ln(2) terms appear here artificially from the absorption of such terms into the factorization and 
renormalization log's, i.e., ln(Q 2 //i 2 ) — ► ln(2Q 2 //i 2 ) — ln(2), see Eq. ()2.35|) . All other expressions 
in Eqs. (j3.43l - l3.47j) remain unchanged, e.g., 

^CS(l) _ rpCS(l) rpCS(2) _ rpCS(2) rpCS(2) _ rpCS(2) /„ ^ 

I j ~ I j ' I F,j ~ I F,j 7 l Gj - 1 G,j ■ 

3.4 NNLO predictions for small and intermediate values of \u\ 

Based on numerical observations on the small and intermediate behaviour of the transition 
form factor predicted by perturbation theory, interesting phenomenological aspects have been 
pointed out in Ref. jHOl- Unfortunately, for these configurations both photons are virtual, cf. 
Eq. and thus the statistics is rather low. Therefore, no measurements has been done yet, 

although they could be possible at the existing e + e~ machines of the Babar, Belle and CLEO 
experiments. In the following we want to add some comments on the pion transition form factor 
in the small and intermediate \u\ regions and to give predictions at NNLO. 

From the representation ()2.34|) it follows that the jth conformal moments for \u\ < 1 are 
suppressed by u 3 . In addition, the hypergeometric functions appearing in the Wilson-coefficients 
are sharply peaked at \u\ = 1 owing to a logarithmical enhancement caused by the ln[(l — 
u>)/(l + u)} term. For fixed |o;| < 1, one finds for growing j an increasing suppression of the 
hypergeometrical functions, in addition to the power-like suppression due to u J . To study this 
behaviour in more detail, we employ the integral representation for the hypergeometrical functions: 



— 2? 1 



2cj 



[l + uY+^i- V 2(j + 2 + | 7j ) 
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1 + LU 



(3.51) 



r(4 + 2j + 7i ; 



dsse 2 e 



|-s 2 +ln(l-e- s2 )-ln(l-o; 2 e- s2 )| 



r(2+j+ 7 ,/2) 2 Jo 

To evaluate this integral for large j, we rely on the saddle point approximation, which is valid as 
long as the condition 



UJ 2 > 1 



(3.52) 



is satisfied. To clearly illustrate the suppression we mentioned above, we write the Wilson- 
coefficients in the form 



d 



Cj = Cj [a a (n),Q/fi, — 



r(j + i)r(j + 2)r(2j + 4 + 7j )^ 



dj \ 2Q 2 



T(j + 2 + 7i /2)r(j + 3 + 7i/2)r(2j + 3) 

Consequently, we have the following normalization for the function Ej(cu = 1|7j) 
()3.38|) . For \u\ < 1, the approximation 



£>|7i)- (3.53) 
1, cf. Eq. 



v^F J2 (j + 1 + | ) VT 



2V1 + VT 



exp j- (j + 1 + ^ In (l + Vl - u 2 ) } (3.54) 



shows an exponential decrease for (j + 1) ln(l + a/1 — a; 2 ) > \/2 . The value of i^(u;|7j) is then 
smaller than 1/2. Note that perturbative corrections due to the anomalous dimensions, which are 
positive and grow logarithmically with j, give a logarithmical enhancement of this behaviour. In 
the case of rather small \u\, the suppression factor is proportional to 2~i~ 1 ^/j + 1, which affects 
even the lowest partial wave j = 2. The suppression is already larger than 80% for given j > 6 as 
long as the inequality 1 — \u\ > 4/(j + l) 2 is satisfied. Increasing will then abruptly increase 
the value of Ej to reach the Ej = 1 limit. To finish this general discussion, we estimate the partial 
wave that will be suppressed by a factor £ = Ej(u) < 1/e depending on u: 



J + l 



-W_i 



4g 2 (l+y / T^I?) ln(l+y / I^ 2 
TTvi— iu 2 



21n(l + vT^ 



- In (4£ 2 /tt) 

4(1 -H) 



(3.55) 



where x) is the product log function that is real valued in the region < x < 1/e. 

To estimate the contribution of the first few non-vanishing partial waves, we first consider the 
conformal moments (ii(P)\Ojj(fif)\Qy ed = QNjBj. The distribution amplitude vanishes at the 
end-points pQ and from this behaviour it follows that NjB~ vanishes at j — > oo: 



IX 



1 — xY for x 



6NjBj ~ j 6 for j — ► oo with e > 0. 



(3.56) 



We want to add that different non-perturbative estimates, based on a lattice calculation, sum 
rules or a model calculation give quite different values of 6N 2 B 2 at a scale Q < 1 GeV, varying 
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from ~ — 1 to ~ +1. Here the lower bound stems from a preliminary lattice calculation [64], while 
the upper one arises from sum rule estimates EE] and is also compatible with previous lattice 
calculations (see jS] for references). There are other estimates that favour a rather small value 
of B 2 . This suggests that the absolute size of the lowest few conformal moments 6NjBj are of the 
order one or even smaller. In the following estimates we consider them of the order one, which 
serves us as an upper bound for the contribution of the jth partial wave to the transition form 
factor. 

In the small \u\ region, i.e., \u\ < 0.4, the lowest partial wave contributes essentially. In LO 
the relative contribution of the second and the forth partial wave with respect to the first one is 
for \u\ = 0.2(0.4) about 0.08(2.3)% and 0.004(0.05)%, respectively. The u 2 term of the zeroth 
partial wave varies in the same order as the relative correction to the second partial wave, which 
is in addition suppressed by a relative factor of 2/3. Thus, in the small \u\ region perturbative 
QCD provides us an (almost) parameter free, factorization scheme independent prediction: 

F 7w (a;,Q)~^CoK(/i r ),g//i r )^l + y + 0(a; 4 )^ for M<0.4. (3.57) 

The phenomenological consequences are obvious, since this prediction is practically independent 
of uj and its logarithmical Q 2 dependence is governed only by the running of the coupling. 

For intermediate values of \uj\, defined as 0.4 < \oj\ < 0.8, the second partial wave contributes 
between 2% and 13%, while the fourth one is at least more than five times suppressed with respect 
to the second one. Increasing |u;| to the value 0.95, the relative contribution of the second and 
fourth partial wave grows to 25% and 10%, respectively, while the sixth (eighth) partial wave 
contributes at the 4%(2%) level. It is illustrative to compare these numbers with the suppression 
arising in the limit \oj\ — > 1 in which the contribution amounts to 39%, 24%, 18%, 14% for 
j = 2, 4, 6, 8 partial waves, respectively. 

As we have realized,only the first two non-vanishing partial waves are essential for an in- 
termediate value of |co>|. It could, therefore, be justified to employ the Taylor expansion of the 
hypergeometric functions at uj = and hence, the transition form factor reads in the CS scheme 

72 



x6N 2 B 2 (fi 2 ) + 0(J 



for 0.4 < \u\ < 0.8, (3.58) 

where the Wilson-coefficient c 2 (a s (/i), ^, is defined by Eq. (J2.35J) to NNLO accuracy and a 
consequent expansion in a s should be done to this accuracy. For \uj\ = 0.8, the higher-order terms 
in 0(uj 6 ) contribute at the 2% and 20% level for the j = and j = 2 partial wave, respectively. 
These contributions can be reduced by a factor of two [four] by taking the order uj 6 [u 8 ] corrections 
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into account. For larger values of \u\, the convergence of the Taylor expansion at uj = is rather 
slow for higher partial waves. For instance, to approximate the third non-vanishing partial wave 
at \uj\ = 0.9 to an accuracy of better than 10%, one has to take into account the first ten non- 
vanishing terms, i.e., up to 0(u; 24 ). 

For ultra-large values of \uj\, let us say \uj\ > 0.95, partial waves with higher conformal spin start 
to contribute with increasing |u;|. However, as we have discussed, as long as we do not reach the 

— > 1 limit, there will be an exponential suppression for higher values of j. Note that this limit 
can never be reached in any experiment at an e + e~ machine, where the mean value of the virtuality 
of the untagged photon is set by the electron mass and there are further kinematical restrictions 
arising from the detector geometry and kinematical cuts. Just for illustration, we would like to 
mention that for |a>| = 0.99 and \u\ = 0.999, the contributions of the 12th and 38th partial waves 
are reduced by a factor 1/e ~ 0.37 compared to their contributions in the limit \u\ — > 1, while 
higher ones start to be exponentially suppressed, since the (j + 1) ln(l + VT — u 2 ) > \/2 condition 
is fulfilled. 

We now present the general result of the photon-to-pion transition form factor for \u\ < 1 in 
its expanded form to the NNLO approximation. Notation analogous to that in Eqs. ()3.42j) will 
be used and the Taylor expansions in uj 2 for the contributing terms of the first five non-vanishing 
partial waves are listed in Appendix [01 At leading order the hypergeometrical functions can be 
expressed in terms of elementary arctanh functions 5 , e.g., 

2 . arctanh(u;) 



1 - 1-u/ 



UJ 



(3.59) 



2 r 2 4 3arctanh(o;) 



15- 13a/ - (5 - Quj 2 + cj 4 )- 



OJ 

and the expansion in uj is given in (jD.2|) . The radiative corrections for uj ^ depend on the 
factorization scheme even for the lowest partial waves. In comparison with Eqs. ()3.43j) and ()3.44j) . 
here we will not include ln(2) terms in the factorization and renormalization log's, i.e. instead of 
2Q 2 we rather employ the underlying scale Q 2 ~ 1/z 2 of the OPE of two currents at the distance 
z. Thus, we have for the MS scheme the NLO result ()3.16|) 

Tf\u,Q/ii f ) = C F (r^{uj) + Tf\uj)\n v^j , (3.60) 

and analogous one for the CS scheme. Note that in comparison to this notation the definitions 
of conformal moments in the |a;| — > 1 limit, Eqs. ()3.32|) and ()3.45|) . differ by a ln(2) proportional 



term: 



lim Tghu) = Tp] + ln(2) t\ 0) v 3 , lim T^ (1) (uj) = T^ (1) + ln(2) T^v, . (3.61) 



3 The result can be expressed in terms of In functions by means of arctanh(w) = 1/2 ln[(l + uj)/(1 — uj)] . 
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Expanding the T^(lu)LN 1 (uj) terms in Eq. ()3.18|) provides after convolution with the correspond- 
ing kernels the desired result in the MS scheme (ID. 8(1 . In the considered order of a; 2 , it coincides 
with the result of Ref. jHOj- The result in the CS scheme can be easily derived by expanding Eqs. 
(j2.18j) in order a s \ 



(3.62) 



where c/ 1 ) and Vj are given by Eqs. (j3.41j) and (|3.34j) . respectively, and s * (uj) is defined by the 
expansion 



( j+l+e,j- 


h2+e 


2lu \ 


\ 2(j+2H 


-*) 


l+u> J 



(1 + 



1 



uo)e 2 + 0(e 3 



2^1 



3+1,3+2 
2(3+2) 



2a; 



1 + W)^ 1 



(3.63) 



The corresponding expansions of and sf in u; 2 are given for j = 2, 4, 6, 8 (because of current 
conservation Vq — and, consequently, j = term does not contribute) given in Eqs. (|P.5)1 and 
(jD.7|) . respectively. The NNLO contribution in the CS scheme we write analogously to Eq. ()3.44j) 

as 



.(2) 



Cj? < CpT, 



CS(2) 



F 1 F,j 



A) T (2)/ N 



x 



C!f ( «fy + 4% - ^» + ^ln I ^ I I - 



,2„(1) 



o 2 



A) 



CS,S/ 



a; 



(3.64) 



- C 



0i 



Q 2 



-,CS(1) 



Q 2 



In comparison to the definitions in Eqs. (|3.46j) . (|3.47|) . and (j3.33J) . one has again to take care of 
ln(2) terms by means of the following correspondences: 



lim^fV) 
limT G C f ( W ) 



,(1) 



lim To „• l a; J 



rg (2) + in(2) rf ] 

Tg (2) + ln(2) if > Voj 
(i) 



V FJ + C) >Vj 



Va s 



l) + ^n(2)) 



(3.65) 



2$ + ln(2) 



if^j' 2 - T^ 1 ' - ~ ln(2) T/> 



nCS(l) 



1 

2 



(o), 

3 "3 



and similarly for the quantities in the CS scheme. Here we have also employed the identity (jP.71 



For Tpf 2 \oj) and T, 



CS(2) 
G,j 1 



o>), the expansion of the Wilson-coefficients ()2.18|) gives 



T°f\u) 



T 



(o), 



.(2) 



,(1) 



W ) \ C FJ ~ \ C j V 3 + V F,j S) \UJ) + -fs) '{U 



^' (2), 



CS(2) 



(a;)=Ti°V){c^-^ a f(a;)}. 



(3.66) 
(3.67) 



The conformal moments proportional to /3q, are obtained from the MS result given by Eq. ()3.20j) . 
The expansion in uo 2 of the term Tgj is given by (jP.9j) . As in the 
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1 case (see the discussion 



of expression (|3.35Jl ). we define the conformal moments of T^(uj) £g> [vp], by Tj (u) VpAoS). Note 
that the "effective" conformal moment i;^ • now depends on u and its expansion in uj 2 is given in 
(|D.10p . Analogously to Eq. ()3.48|) . Vg-iui) is provided by 



T r ra(1) y-rf»y 



V M = W w ) + ^(o), : — > ( 3 - 68 ) 



and the expansion can easily be obtained by means of Eqs. ()3.62|) . (jD.5|) . and (|D.8|) . 

The fa proportional NNLO terms in the CS scheme are obtained by making the replacements 

Tgj H - T^ (2) ( W ) = if ( W ) {eg - %j sf ( W ) } and - v g = %j if ( w ) , (3.69) 

while the other terms remain the same as in the CS scheme. 

4 Radiative corrections to the photon-to-pion transition 
form factor 

This section is devoted to a model independent study of radiative corrections to the pion-to- 
photon transition form factor in the case of one quasi-real photon — > 1) and in the small 
and intermediate \u\ -regions. We also illustrate how the perturbative QCD approach to exclusive 
processes can be tested in a novel way by a sum rule and how the two lowest non-trivial conformal 
moments of the pion-distribution amplitude could be extracted from experimental data in the 
intermediate |o;|-region. 

In Section B~T1 we briefly review the features of the radiative corrections to the first few confor- 
mal moments of the hard-scattering amplitude in the MS and CS schemes to NLO. We point out 
that asymptotic formulas with respect to the conformal spin j + 1 provide a very good approxi- 
mation of these moments in question for a rather low value of j > 4. As a byproduct, we propose 
a simple method for reconstructing the amplitude from its conformal moments, which is outlined 
in Appendix |Ej 

In Section \A. 21 we present the numerical values of the NNLO corrections to the first five non- 
vanishing conformal moments of the hard-scattering amplitude in the CS and CS schemes. We 
point out their general features and discuss different possibilities for treating the (3q terms. In 
particular, we consider the lowest conformal partial wave and compare its contribution to the 
photon-to-pion transition form factor with experimental data. We study the influence of radiative 
corrections to the sum rule and show that higher-order corrections will interfere only slightly in 
the extraction of the two lowest non-trivial conformal moments of the distribution amplitude. 
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Table 1: First moments of Tpj, Tpj, and T^J for |a>| = 1 with /A = 2Q 2 and for \u\ = 0.8 with 

n- Q 2 . 





j 





2 


4 


6 


8 


10 


12 


14 


16 


18 


\u\ = 1 
li) = 2Q 2 


T (o) 

T (l) 

F,j 

rpCS(l) 

J F,j 


1.5 
-3.75 
-2.25 


0.58 
1.2 
1.91 


0.37 
2.14 
2.52 


0.27 
2.38 
2.58 


0.21 
2.42 
2.51 


0.17 
2.39 
2.41 


0.15 
2.33 
2.3 


0.13 
2.26 
2.19 


0.11 
2.19 
2.09 


0.10 
2.12 
2 


\u\ = 0.8 

v} = Q 2 


T (0) 

T (l) 

F,j 

rpCS(X) 

2 F,j 


1.19 
-2.02 
-1.78 


0.15 
0.14 
0.16 


0.03 
0.1 
0.1 


0.01 
0.03 
0.03 


0.01 
0.01 













4.1 Features of radiative corrections at NLO 

Let us first compare the NLO corrections in the MS and CS schemes for \u>\ — > 1. The first ten 
non- vanishing moments Tp! and T F j , given by Eqs. ()3.32|) and ()3.45|) . are shown in Tabled 
From this table we realize that the main difference between the two schemes is in the first two 
moments which differ by about 50%, and in a seemingly faster decrease of the moments in the CS 
scheme for large j. Indeed, in the large j asymptotics the leading terms are 

2j + 3 f c/i a 9 



/, - / (j + 1)0' + 2) 



tCS(i) _ 2j + 3 

a +2) 



t fT r.w/, ^ wi+j 



{2S 1 (l+ J y-- + 0((3 + 2r 1 )Y (4.1) 

C(2)-31n(2) + 0((j + 2)- 1 )}. 



Si(l+j) + ^+4m(2) 



9 
2 



Taking into account the large j asymptotics of the Si functions, given by 

limSi(l+j) = ln(2+j)+7E, (4.2) 

the ratio t£?^ I^fI slowly approaches 1/2. The difference is caused by the (infinite) resummation 
of off-diagonal terms in the MS scheme. The asymptotic formulas ()4.1)1 have a relative error of less 
than 2% already for j > 4. Thus, by knowing a few lowest moments and their asymptotics we gain 
a complete insight into the radiative corrections for \u\ — 1. In Appendix |E] we use this result to 
make an approximate reconstruction of the hard-scattering amplitude in the momentum fraction 
representation from its conformal moments. The consequence of the logarithmical behaviour in 
Eq. ()4.2j) is obviously an increase of radiative corrections with growing conformal spin. It is shown 
in Table H3 that already for j = 8, radiative corrections are of the size of 80% for a s /7r ~ 0.1 (i.e., 
[L r fa 2 GeV for one-loop a s with rif = 3). From this point of view, one might conclude that 
perturbation theory breaks down for rather large values of j. Fortunately, higher conformal spin 
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Table 2: The first five non- vanishing Wilson-coefficients appearing in the perturbative expansions 
of Tp s and TP S with respect to a s /rr at NNLO accuracy for \u\ = 1. The results are obtained 
employing fij = jx 2 = 2Q 2 . 
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13.56 


11.33 
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0.27 


1.72 


1.25 


-0.58 


15.17 


12.06 
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0.21 


1.67 


2.54 


-0.55 


15.68 


12.05 



contributions are suppressed by the non-perturbative input, see Eq. (J3.56J) . and so perturbative 
QCD remains applicable. In the photon-to-pion transition form factor there might be also a 
cancellation of the lowest partial wave with the remainder, which is due to their relative minus sign. 
Of course, the net contribution of radiative corrections depends on the model of the distribution 
amplitude itself. 

With decreasing |a>|, higher conformal partial waves are starting to be exponentially suppressed, 
and as we have shown in Section l3~4l radiative corrections logarithmically enhance this suppression. 
Also note that off-diagonal contributions to each partial wave, which are relatively suppressed by 
powers of uo 2 with respect to the diagonal ones, are becoming small. If we approach the equal 
virtuality case, i.e., u = 0, only a factorization scheme independent constant, arising from the 
lowest partial wave, will survive. Thus, by decreasing \u\ the differences between the MS and 
CS schemes must be washed out. In Table [^we illustrate these effects for = 0.8. For the 
two lowest non-vanishing partial waves the differences between these two schemes is reduced to 
about ±14% and for higher ones below 2%. Also in the CS scheme the contributions from the 
functions s®(u)) = 0(u 2 ), cf. Eq. (j3.63J) . are power suppressed. So one expects from Eq. ()3.62|) 
that radiative corrections due to the normalization factors are the essential ones, but with one 
exception. Since the coefficient is relatively small compared to the anomalous dimension ^\ 
0{uj 2 ) corrections remain important for the second partial wave in the intermediate \u\ region. 



33 



4.2 Predictions to NNLO accuracy 
4.2.1 The quasi-real photon limit 

We now turn to the discussion of NNLO effects in the CS and CS schemes starting with the limit 

CS(2) ^08(2) 



\U\ 



r CS(2) 



1. In Table El we present the numerical values of the Wilson-coefficients T F , T ( 



TP], and T, 



CS(2) 



■.a •, auu corresponding to Eqs. ()3.46|) . (|3.47|) . (|3.33|) . and ()3.49|) . respectively. 



The values of Tp^ 2 \ 



Tq S ^ 2 \ and Tg^ (Z> have been obtained by means of the NNLO result for the 



-,CS(2) 



deep inelastic scattering structure function gi [43j. 

Let us investigate in more detail the contribution of the lowest partial wave to the transition 
form factor, which is scheme dependent for u ^ and for \u\ — 1 reads: 



in the MS scheme: 
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For a s (nl)/ir = 0.1, the ratio of the NLO to the LO contribution is —17% in the MS scheme, and 
—10% in the CS scheme. This difference arises from the fact that in the MS scheme off-diagonal 
terms of the hard-scattering amplitude are resummed. In Eq. ()4.3|) we see that the ln(2Q 2 /^ 2 ) 
term is rather small compared to the ln(2Q 2 //U 2 ) one. This is even more the case in the CS scheme, 
while in the CS scheme the ln(2Q 2 //i 2 ) term vanishes completely, since all off-diagonal entries in 
the NLO evolution have been removed. The sign alternating series of the (3q non-proportional 
terms is due to the Sudakov effect, see Ref. (2H| for a detailed discussion. 

Since factorization scale changing effects in the hard scattering amplitude are for the lowest par- 
tial wave quite small and since they will be compensated by the evolution of the non-perturbative 
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part (see Section [23} , we set /j, 2 = 2Q 2 in the following and discuss the scale setting of the residual 
ji r dependence. First, let us equate fif, = 2Q 2 : 



F'fTT (Q) 



y/2f„ 
2Q 2 



a s (2Q 2 ) [7.23] a 2 s (2Q 



71 L 5.14 J TT 



2\ 



r'2 



0{al) 



for {Sq -scheme. (4.6) 



Hence, for a s (fi 2 .)/7T = 0.1, the ratio of the NNLO to the LO contribution is —7.2% and —5.1%, 



and the ratio of the NNLO to the NLO contribution (the measure of the convergence of the 
perturbative QCD expansion) is ~ 70% and « 50%, in the CS and CS schemes, respectively. 

The main part of these rather large NNLO contributions arises from the /3o-proportional term. 
Owing to the off-diagonal parts, it is larger by about a factor of two in the CS scheme than in the 
CS scheme. It is appealing to resum this large contribution by the Brodsky-Lepage-McKenzie 
proposal [HZj (for application to exclusive processes see also jHHl); hi which all terms proportional 
to (3 are absorbed in the coupling by the scale setting /i r = ^slm- 



VjjA f a s (/i B LM) nno «s(/v) , /v :' M 

2Q 2 Y ix +U ' y2 n 2 



Fm{Q) = ^\l- SK T J +0.92^^ + 0(ai)} , (1.7) 

with 



2 ^ 2 f 1/37.431 r fCS , , A , 

^BLM = 2g 2 |^ 14 78 | for [—-scheme. (4.8) 

The ratio of the NNLO to the NLO coefficient is now only minus one and reflects the Sudakov 
effect in the conformal theory. However, as we realize, combining the COPE result with the MS 
result of the /^-proportional piece induces a rather low scale. For instance, for 2Q 2 = 4 GeV 2 we 
have /^blm ~ 0-1 GeV 2 in the CS scheme and hence the non-perturbative behaviour of the coupling 
is needed. If we completely remove the off-diagonal terms, the BLM scale squared is enlarged by 
a factor of 2.7 and is now closer to that in the MS scheme |XQJ given in Table |U What one is 
actually doing here is to combine perturbative QCD with speculations about the non-perturbative 
behaviour of the QCD coupling and so, strictly speaking, one is leaving the perturbative ground 
on which the whole analysis was based. However, one advantage of this proposal is that the 
result predicted by conformal symmetry is recovered if we consequently assume a hypothetical 
fixed-point of the /^-function during our considerations. What we in fact do by the freezing of the 
coupling, is to assume that this non-perturbative fixed-point is at Q 2 = 0. 

In Fig. ^ we compare the experimental data from the CLEO experiment with the prediction 
arising from the lowest conformal moment, the only one that survives in the asymptotic limit 
Q 2 — > oo. Without further considerations, we assume, as has also been done in the employed 
method of data extraction (see Section 3 of Ref. [H]), that the quasi-real photon limit has been 
reached. The prediction for the asymptotically large Q 2 in this |a;| — > 1 limit is displayed as 
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Figure 1: The contribution of the first partial wave to the scaled photon-pion transition form 
factor 2Q 2 F^ n (u = ±1, Q) is shown in LO (dotted), NLO (dash-dotted), and NNLO (solid and 
dashed) accuracy for the CS scheme. The renormalization scale has been set to /i 2 = 2Q 2 (solid) 
and to the BLM scale /i 2 = 2Q 2 /14.78 (dashed). The data are taken from Refs. [T^ITH]. 

a dotted line. As we have discussed, radiative corrections reduce the size of this prediction for 
realistic values of Q 2 . The dash-dotted line represent the NLO and the solid line the NNLO 
prediction for the standard MS definition of the coupling in the corresponding approximation 
with the normalization a s (fi r = Mz) = 0.118. The dashed line is the prediction for the BLM scale 
setting and accumulates non-perturbative effects by freezing the coupling to a s (fi r = 0) = 0.6 by 
adding an effective gluon mass. Let us note that the scale setting ambiguities at 2Q 2 = 1 GeV 2 
are of the order of 20% if we vary the renormalization scale /z 2 from 0.5 GeV 2 to 2 GeV 2 . This 
ambiguity can be further reduced by going to higher orders in a s . 

This comparison shows that theoretical uncertainties due to higher order radiative corrections 
and scale setting ambiguities are much smaller than the error of the experimental data. At 
larger value of Q 2 there is no significant contribution of higher partial waves. However, there is a 
significant discrepancy of the results in Fig. 1 in the region 0.5 GeV 2 < 2Q 2 < 2 GeV 2 , which may 
indicate the presence of higher partial waves. Since evolution effects in this kinematical region are 
rather strong, we could employ them to pin down the size of higher partial waves. However, the 
question arises: Can we in this kinematical window rely on the leading twist result? 

By considering the size of power suppressed contributions, we will now argue that the answer to 
this question is positive. Since the (local) matrix elements of any operator appearing in the OPE 
can only be built with the momentum four- vectors P, Lorentz covariance immediately tells us that 
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Table 3: The ratio of NLO to LO and NNLO to NLO radiative corrections in units of a s /iT for 
\u\ = 1 and y? T = \A — 2Q 2 in the MS, CS, and CS schemes. 
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2T (0) 
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9.21 


7.39 
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7.64 




7.93 


10.56 
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power suppressed contributions are of even power in Q. The only dimensional parameters that 
can appear are the mass of the pion ~ 0.14 GeV and the QCD scale parameter A ~ 0.2 GeV. 
Assuming that multi-partonic correlation functions will not have a strong numerical enhancement, 
we expect that the contributions proportional to m 2 will provide a relative correction of the order 
of 0.02 GeV 2 /2Q 2 . The size of the remaining non-perturbative corrections arises from ambiguities 
in summing the perturbative series and can be estimated in the framework of renormalons. In 
the conformal scheme, we might again borrow the results from the analysis of the coefficient of 
the structure function g%, which gives for the lowest moment an uncertainty of a similar size, see 
Rcf. [69 1 and references therein. Thus, we might conclude that the non-perturbative effects to the 
contribution of the zero partial wave are smaller than 10% at a scale of 2Q 2 ~ 1 GeV 2 . Certainly, 
this rather optimistic speculation should be confronted with other methods used. The lowest 
conformal moments of twist-four quark-gluon-quark operators have been obtained by means of 
QCD sum rules [7UJ d]. Here it has been found that, relative to m 2 , a certain matrix element is 
numerically enhanced by a factor of ten. This certainly would strike our point of view and indicate 
that the so-called Wandzura-Wilczek approximation, in which higher multi-parton correlations are 
neglected, fails. Renormalon induced corrections have been studied in a model-dependent way in 
Ref. [Hj, where their relative size has been estimated to be 0.2 GeV 2 /2Q 2 . Note, however, that in 
this analysis excitations of higher conformal partial waves have been included. Making it short, 
we stress that the estimate of higher twist-contributions has to be made in a consistent framework 
that is set by the scheme in which one started. Combining estimates from different approaches is 
a popular but rather awkward procedure. 

In this process with a quasi- real photon higher partial waves are summed. Even if the |a;| — > 1 
limit is not reached in the experiment, a rather large number of terms will contribute. Without 
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Table 4: Ratio oblm = /^blm/^Q 2 of the BLM scale squared to 2Q 2 and the ratio of the NNLO 
to the NLO coefficient in units of a s /n for = Hblmi A 4 / = 2Q 2 , and |u;| = 1. 
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any knowledge about the shape of the distribution amplitude, it is a rather vague assumption to 
truncate this series by hand to extract the values of the lowest partial waves from the normalization 
of the pion-to-pion transition form factor. Figure ^ clearly shows that the dominant contribution, 
at least for 2Q 2 > 2 GeV 2 , arises from the lowest partial wave and the remainder is small. The 
fact that the contributions of higher partial waves cancel each other is not excluded and it remains 
a claim that the asymptotic shape of the distribution amplitude is established by experimental 
data. In principle, one can gain more information on the remainder if one also employs the 
evolution of the distribution amplitude. However, even if rather high precision data are available, 
the deconvolution problem is not easy to solve. As we have already mentioned, at NLO the 
perturbative correction will increase with growing conformal spin. The same tendency can be 
read off from Table H3 also in NNLO, where the fio proportional term is the dominant one. This 
is also reflected by the decrease of the BLM scale as shown in Table HJ where we can also see 
that the remaining corrections at NNLO are moderate. Note that the BLM scale is rather low for 
2 < j in the MS and CS schemes, which is due to off-diagonal terms. 

4.2.2 What can we learn from the small and intermediate \u\ regions ? 

As it has been clearly spelt out in Ref. [3U] and explained in a more general way in Section 13.41 
the small |u;| region is suitable for a novel test of the perturbative QCD approach to the class 
of exclusive light-cone dominated processes. As we noted in Section [4.11 for decreasing \u\ the 
differences between different schemes will decrease too. This is illustrated for \u\ = 0.8 in Tables 
El and IH1 For the lowest [second] partial wave we have about a —40% [10%] effect at NNLO 
compared to NLO or a —4% [< 1%] correction compared to LO for fj, 2 = ft 2 = Q 2 and a s /n ~ 0.1. 
Altogether, we find ~ —15% reduction of the LO prediction for the lowest partial wave and an 
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Table 5: Same as in Table El for \u\ = 0.8 and = fi 2 — Q 2 . 
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Table 6: Analogous to Table |U for the ratio oblm = /^blm/Q 2 ' where fj/j = Q 2 and \u\ = 0.8. 
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increase of about 7% and 33% for the second and fourth one, respectively. The main part of the 
NNLO corrections arises from the (3q proportional term. Its absorption in the running coupling via 
the BLM scale setting prescription again requires the knowledge of the non-perturbative behaviour 
of a s . Table El shows that then a sign change occurs at NNLO, where the BLM scale for the second 
partial wave is quite low and its remaining NNLO correction is rather large. 

Certainly, the resummation of the (3 proportional corrections is associated with a new input 
that is not well known. Thus, in the following discussion concerning the extraction of non- 
perturbative conformal moments of the distribution amplitude we prefer the naive scale setting 
prescription p? r = ji 2 = Q 2 . In panels (a) and (b) of Fig. [2] we display the cu-dependence for 
the scaled photon-to-pion transition form factor evaluated in the CS scheme at LO and NNLO, 
respectively. One clearly sees that the prediction is almost independent of uj for |a?| < 0.2 and 
only a negligible dependence arises for 0.2 < \u\ < 0.4. Radiative corrections will only shift this 
prediction downward. Note that this shift will slightly increase, if we go to higher orders of a s . For 
the lowest partial wave they can be taken from the calculation of the radiative corrections to the 
Bjorken sum rule, which are evaluated in the third-loop approximation [72] and roughly estimated 
at four loops [TS]. Consequently, confronting these predictions with experimental measurements 



39 



. 18r 





0.1 



0.2 



0.6 



. 8 



0.4 0.6 0.8 0.2 0.4 

(a) ^ (b) 

Figure 2: The ^-dependence of the scaled photon-to-pion transition form factor 2Q 2 F 17T (uj, Q) at 
Q 2 = 2 GeV 2 in LO (a) and in CS at NNLO (b) for three different values of B 2 = {0, -0.5, 0.5} is 
shown as solid, dashed and dash-dotted lines, respectively. The spread of the corresponding lines 
display the sensitivity of the predictions with respect to the parameter £> 4 , which is equated to 
—0.25 and 0.25 for the lower and upper curve, respectively. 



would provide either a novel test of perturbative QCD or an insight into the size of higher-twist 
contributions. To enhance statistics, one can even integrate over the small \u\ region: 
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(4.9) 



(~ 200)^1 + O(al) + O (m 2 JQ 2 , A 2 /Q 2 ) \, 



where u; cut < 0.4 and rif = 3. 

If we vary the unknown B 2 and B$ parameters in the range that is suggested by non-perturbative 
estimates, i.e., —1/2 < B 2 < 1/2 and —1/4 < B 4 < 1/4, it can be realized that in the interme- 
diate \u\ region the perturbative QCD predictions start to be dependent on the value of B 2 and 
for larger |u;| even on that of B±, while higher partial waves can be safely neglected. Here it 
is important that radiative corrections do not spoil a possible extraction. They rather shift the 
curves in the whole region and slightly enhance the spread of the curves associated with B4. This 
is caused by the fact that the radiative corrections to the lowest two partial waves are negative and 
approximately of the same size, while they are positive in all other cases. Since the gap between 
the curves for both extreme cases, where B 2 = ±0.5 is about 30% or even larger, we might expect 
that these curves can be distinguished by a measurement. This would also test the reliability 
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of non-perturbative methods. Certainly, a precise extraction of B 2 or even B 4 also requires high 
precision measurements with high statistics. Assuming that such data are available, the extraction 
of the non-perturbative parameter itself is rather simple. First, a measurement in the small \u\ 
region should be confronted with the parameter free prediction that depends only on the running 
of the coupling. As we argued, we do not expect higher twist corrections to be essential. If that 
were experimentally established, a simple two parameter dependent fit 



F y7r (io, Q) — f n 



T (lu, Q) + yT 2 (^, Q)B 2 (Q) + g T 4 (u, Q)B 4 (Q) 



(4.10) 



could be used to extract the two conformal moments B 2 (Q) and B 4 (Q) as long as the statistics is 
high enough. Moreover, a consistency check is provided by the Q-dependence of these parameters. 
Note that in the CS scheme the mixing between different partial waves is caused by evolution at 
NNLO, while in the MS-scheme the mixing appears already at NLO in both the hard-scattering 
amplitude and the evolution of Bi. In the CS scheme the mixing appears first in the evolution to 
NLO accuracy. 

Finally, we want to comment on the size of evolution effects which are caused by the off- 
diagonal entries in the anomalous dimension matrix. We numerically observe that the conformal 
symmetry breaking terms at NLO, (compare the MS scheme results with the CS ones in Tabled 
as well as at NNLO, see Table EI) can provide an enhancement of the corresponding corrections up 
to 50%. A similar relative effect one would naively expect from the evolution due to the unknown 
mixing in the CS scheme arising at NNLO. However, since it does not appear at the input scale 
Qq due to the initial condition (j2.42j) . this mixing effect is in fact small j3U]- For instance, at NLO 
in the MS and CS schemes this mixing effect in the lowest partial wave goes up to -1.3% and 2.3%, 
respectively, for the evolution from Q = 0.5 GeV to Q = 20 GeV and \uj\ = 1. Note, however, 
that in the MS scheme cancellation appears in the off-diagonal terms between the C\ and CfPo 
proportional parts and that the contribution from only the C\ proportional term might be of the 
order of 2%. This number should be compared with the correction in the hard-scattering induced 
by the off-diagonal terms, which is of about 7%. As we discussed above, these corrections are 
reduced for \u\ < 1. Since the ratio of diagonal entries in NNLO to those in NLO is smaller than 
3a s /n for the first five even diagonal terms, we might argue that the ratio of off-diagonal entries in 
the MS scheme is of the same size. Assuming so, one would expect that unknown mixing effects, 
induced by A jfc due to evolution are reduced by a factor 1/3 or even more with respect to those 
at NLO. Thus, we expect from the numbers given above at NLO that the contribution of Aj k 
is smaller than 0.5% in the MS scheme. In the CS scheme, the contribution of the log term in Eq. 
(J2.53|) can be estimated by comparison with the NLO effect in the CS scheme. Here one finds a 
relative contribution smaller than a s (Q) ln[Q/Qo]/7i. For the example discussed above that would 
produce a mixing effect smaller than 1%. Therefore, for higher values of the input scale, e.g., 
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Qo ~ 1 GeV, we expect a rather tiny mixing in the CS scheme. 

5 Conclusions 

We have employed conformal symmetry in the perturbative sector to evaluate the NNLO cor- 
rections for the pion production through two-photon fusion. The requirement of a manifestly 
conformal invariant result partly removes the ambiguities arising from the factorization. How- 
ever, the ambiguities are retained in the scheme dependence of the forward Wilson-coefficients 
and anomalous dimensions and in the treatment of the conformal symmetry breaking induced by 
the trace anomaly, proportional to the /3-function. The latter ambiguity has been studied here 
in two alternative schemes: (i) combining the conformal predictions with the MS result and (ii) 
improving the partial wave decomposition of the conformal invariant theory by the renormaliza- 
tion group equation. The second possibility minimizes the mixing of partial waves and gives us 
an almost good quantum number, namely, the conformal spin. For decreasing values of \u\, the 
differences between these schemes are removed, since off-diagonal terms are suppressed by powers 
in u 2 . 

As it has been known before, for = 1 NLO corrections can be considered to be small for 
only the two lowest partial waves, since the NLO corrections logarithmically increase with the 
conformal spin. This behaviour is analogous to the large j behaviour of the Wilson-coefficients in 
DIS arising from soft gluon configurations. The effect is manifested by ln(l — x) - terms that are 
associated with factorization log's and, consequently, are absent in the lowest partial wave in the 
conformal schemes. Other ln(l — x) terms are related to the Sudakov effect and are manifested 
in a sign alternating series for (3 = 0. However, the numerical study has shown that the NNLO 
corrections are dominated by the /3q proportional term, as expected. In general, this term is rather 
large compared with the NLO coefficient and thus the BLM scale is rather low, which drives the 
coupling in the non-perturbative region. 

We have compared the NNLO predictions with the existing data for the quasi-real photon 
case and have found that for 2Q 2 > 2GeV 2 the contribution from the lowest partial wave is 
compatible with the data. The deviation below this scale is induced by non-perturbative effects 
or by the contribution of higher partial waves and it requires a deeper insight in power suppressed 
contributions. Although, there is no doubt in the literature that the CLEO measurement can be 
analyzed in this limit, we should state here that partial waves with sufficiently large conformal 
spin will be exponentially suppressed. This can affect the analysis only if the matrix elements of 
conformal operators with rather large conformal spin will contribute. On the other hand, we know 
that the net effect of all partial waves with the conformal spin j + 1 > 3 is small. Obviously, this 
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does not necessarily mean that the matrix elements themselves are small and therefore we cannot 
say that the asymptotic form of the pion distribution amplitude is experimentally established. 
Also, having a look at non-perturbative results from other methods given in the literature, a 
strong statement that the asymptotic form is suggested by those estimates cannot in fact be 
made. 

These problems we have spelled out, can be separately studied apart from the — ► 1 limit. 
Indeed, in the small |u;| region, perturbative QCD predicts a sum rule that has the same status 
as the Bjorken sum rule in deep inelastic scattering, evaluated at order o? s . A first test of this 
sum rule might be possible with existing e + e~ machines and would offer us a first insight into 
the size of power suppressed contributions for exclusive processes from experimental data. We 
expect that such contributions will turn out to be small. If this should be established experimen- 
tally, one might attack the extraction of the first and second lowest conformal moments of the 
distribution amplitude. This is an important task, since it would open a window to testing the 
reliability of non-perturbative methods applied to exclusive quantities. Having in mind that the 
collinear factorization, applied here to the photon-to-pion transition form factor, is also adopted 
for the analysis of exclusive B-physics, it is timely to confront such methods with experimental 
measurements. 

Let us finally give a short outlook for the application of the conformal approach to other 
processes. After a simple replacement of the decay constant and matrix elements the NNLO result 
obtained can be used for the analysis of rj production, i.e., its flavour octet component. Moreover, 
the formalism can be extended in a straightforward manner to the 7/-to-photon transition form 
factor. Guided by the large j + 1 asymptotics of the conformal moments, it is also possible to 
reconstruct the hard-scattering amplitude in DVCS. The reliability of this technique can be tested 
at NLO and partly also in NNLO, i.e., for (3 proportional terms. We also want to add that one can 
go one order further in a s in the approximation of the first few conformal moments of the hard- 
scattering amplitude, since we can borrow the forward Wilson-coefficients from the non-singlet 
sector of the deep-inelastic structure function F 3 , evaluated in N 3 LO 58 . 
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A Structure of the hard-scattering amplitude in the MS 
scheme 



Here we prove that the LN(a>, x) terms, appearing in Eqs. ([3.18< - ET2*Tj) . are related to the factoriza- 
tion log's. We do not distinguish between the renormalization and the factorization scales, since 
the difference appears only in the cu-independent term ln(/x r //x^). The hard-scattering amplitude 
is given by the sum over all Feynman diagram T contributions 



— T(uj,x,Q/ii) = J2Tf(u),x,Q//j,), (A.l) 

W T 



where we have rescaled the individual contributions by Q 2 to have dimensionless amplitudes. Each 
of these contributions is given as a product of propagators and vertices integrated over the virtual 
loop momenta. The two photon vertices are connected by a chain of quark propagators S and 
quark-gluon-quark vertices V: 

#([1 - 2x]P/2 + q- l)Vi - 2x}P/2 + q - I + h) . . .$([1 - 2x]P/2 + q-l + ... + k n ), (A.2) 

where P = q% + q 2 and the large momentum q = (q\ — 92) /2 flows only into this chain. Momentum 
conservation requires that I + I' — Z)"=i ki, where I and I' is the sum of virtual momenta flowing 
into the first and flowing out of the second photon vertex, respectively. Interchanging the two 
photon vertices will give the crossed contributions with \u\ — ► —u. Obviously, there are further 
propagators that depend only on the virtual momenta and xP or (1 — x)P, but not on q. In- 
troducing the Feynman- Schwinger representation for the propagators, integrating over the virtual 
momenta and making use of the on-shell condition P 2 = 0, give us the following representation 
for the regularized contribution: 

f 1 f Zl rzm-i a 2e Tjr(z\e) 

lo Jo Jo Q 2e [xuBi(z) + (1 - x) ujB 2 (z) + 1] + 

Here e is the dimensional regularization parameter and the functions <8j(V) depend on the Feynman- 
Schwinger variables z = {z±, . . . , z m } with n < m. We introduce the new variable 

y=\- X -B 1 {z)- l -^B 2 {z) (A.4) 



and write 



1 fi 2e 

Tj:(u,x,Q/n)= I dy-^-^ ll+e Vjr(y,x\e) , (A.5) 



lo [1 - (2y - l)u 

where the unrenormalized convolution kernel is defined by 

- 1 , C Zl , rZm-l / 1 1 — X 



VAy, x\e) = j o dz x J o dz 2 ... j o dz m 5 [y - - + -B x {z) + — B 2 (z)j T r {z\e) . (A.6) 
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The renormalization procedure provides the factorization log's ln(Q 2 / )i 2 ) , which always appear in 
combination with LN(a>, x) terms. The factorization theorem tells us that after resummation of all 
Feynman diagrams the corresponding convolution kernel is given by the evolution one. Obviously, 
the log independent terms can also be represented as convolution. The support of all these kernels 
is known and follows from the restrictions of the Bi(z_), which are obtained from their definition 
and the topology of Feynman graphs (see, for instance jZHEj). 

B Consistency check with the forward- limit results 

In this section we present a consistency check between the results for the non-singlet coefficient 
function of the DIS polarized structure function g x [13] and the hard-scattering amplitude of the 
pion transition form factor ED]. The former quantity is known to NNLO, while for the latter 
one discussed in Sect. 3.1 the calculation has been performed up to /5 -proportional NNLO terms. 
Both results have been obtained in the MS scheme. Making use of the fact that both quantities, the 
photon-to-pion transition form factor and the structure function g x , are defined by the two-photon 
amplitudes belonging to a general class of the scattering amplitudes for the two-photon process 
at light-like distances, we are able to transform the results for the photon-to-pion transition form 
factor to the results for g 1 . 

The general scattering amplitude for the two-photon processes is given by the time-ordered 
product of two electromagnetic currents sandwiched between the in- and out- hadronic states with 
momenta Pi and P2, respectively. Using the notation q = (qi + ?2)/2 (qi and q2 are incoming and 
outgoing photon momenta), P = Pi + P2, and A = P2 — P x , the following generalized Bjorken 
region can be defined |H |3J5] : 

v = P ■ q — > 00 and Q 2 = —q 2 — > 00 , (B.l) 

with the scaling variables 

u = — and 1] = . (B.2) 

Q z v 

In the forward case, corresponding to DIS, 1/u can be identified with the Bjorken variable xbj 
and 77 vanishes, while for the two- photon production of a hadron rj = 1. The relations between the 
non-forward ER-BL kernels and forward DGLAP kernels were extensively studied and derived in 
Ref. jlj, while in Ref. [HH] the consistency between the transition form factor and g x results has 
been reported up to NLO. Here we explain in more detail the technical side of these consistency 
checks and extend them to /?o-P ro P°rtional NNLO terms. 

The building blocks of the hard-scattering amplitude for the photon-to-pion transition form 
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factor T(u,x,Q, fj,/) given in Eqs. (jl.2|) . ()1.3|) and ()3 .16113. 23*)) can generally be written as 

A^H",*) = i f- 1 n ■ ®*(y-s), 

1 — u;(2y — 1) — ?e 

= — r — ^tt^— +, (B.3) 

1 — c^(2y — 1J — is 

where n — 0, 1, 2. Note that we have reintroduced the term, originally present in the definition 
of Feynman propagators. For the kernels that appear in Eqs. (j3.16H3.23"|) we use a generic symbol 
v(y,x). Furthermore, for the kernels of interest, given in (j3.7H3.8j) and (|3.11|) . the function v is of 
the general form 

v(x, y) = 6(y - x)f(x, y) + { X X ) . (B.4) 

We have to extend our restricted non-forward kinematics to the whole kinematical region. The 
extension of the BL kernels v to the whole x,y region (— oo < x, y < oo) 0] is accompanied by 
the change of the 9 function as 

v ext (x,y) = v(x,y) le{y _ x ^ e{l ^ )e ^ )sisn{y) . (B.5) 

Furthermore, the dependence on rj has to be restored and one performs the following change of 
variables: 

l+t/y l + t'/ V 
x -> , y -> , u^rju. (B.6) 

The definition of the distribution amplitudes, with which the hard-scattering amplitude is convo- 
luted, as well as the definition of its generalized counterpart, introduce the restriction — 1 < t < 1. 
After examining the 9 functions in the kernels (jB.5j) and taking into account \rj\ < 1, one obtains 
— 1 < if < 1. The building blocks of the generalized two-photon scattering amplitude are thus 
obtained and they are of the form 

a , 1 

Ai{u,r],t) = - 



tot — is 



;i ln n (l - cot' - ie) 1 

A 2 , n (u,r ] ,t) = / dif — K - — 

' i 1 — ujv — ie 2rj 



v ' 



(B.7) 



It is easy to see that relations (|B.3)1 indeed represent the r\ = 1 limit of (|B.7|) . 

The forward limit (i.e. the forward Compton amplitude) corresponds to rj — > and, due to 
the optical theorem, the non-singlet coefficient functions of the DIS polarized structure function 
gi (contributing to the total cross-section) are determined by taking the imaginary part of the 
forward amplitude. Taking into account the definition of gi and its coefficient functions (see jSj), 
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one arrives at the following recipe 6 for the building blocks of the non-singlet coefficient functions 



C N8, 



A\ 9l) (z) = -Im 



CO 



71 



lim Ai(u, rj, t) 

?7->Q 



I u=l/z,t=l 



(B.8) 



The 77 — >■ limit of the extended ER-BL kernels [3] results in the corresponding DGLAP kernels 
P of the general form 



P(z) 



P(z) = 0(2)0(1 - z) p(z) with 
It is straightforward to derive 

^tfi + t'/v i + t/v 



P{z) - 5(1 - z) / dz'P(z' 



1 



lim 

ri-*o 2r] 



sign(t) 



P 



with p given by 



t 



P 



lim — 

f7->0 27] 



f 



l+t'/r] l + t/r] s 



f 



i -*7»7 i -tii) 



j + 



The imaginary part of expressions (|B.7J) is obtained by making use of 

1 



Im- 



7r5(l - 



(B.9) 



(B.10) 



(B.11) 



fB.12) 



1 — t'u; — 2£ 

and for more complicated functions, containing ln n (l — t'u — ie) (n = 1, 2), we derive the following 
decompositions 

ln(l - cut' -is) 1 1 



1 — cut' — ie 



ln 2 (l -out' -is) 
1 — cut' — ie 



1 — s)+ 1 — st'uj — ie 
/ln(l — s) — Ins 



1 - s 



<8) 



1 — st'cu — ie 



(B.13) 



Alternatively, for the imaginary parts of expressions, containing logarithms, one can refer to [TS] . 
Finally, we present the results relevant to the NNLO calculation: 



A[ 9 ^ S \z) 

4%^\z) 

4*r s) w 



4^' NS) ( 



6(1 -z) 







Z' J \z' - Z) 



P(z') 



- 1-- 2 



ln(^' — z) — In z 



P(z') 



(B.14) 



V z' — z ) 

6 Factor l/(2w) comes from the dispersion relation, and the additional factor of 2 from the definition of g±. The 
origin of the factor lu in Eq. (|B.8|I lies in the fact that the transition form factor is scaled by Q 2 , while the forward 
Compton amplitude is scaled by P ■ q. 
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Hence, the building blocks for the hard-scattering amplitude of the photon-to-pion transition 
form factor, A^ 1 (u,x) given in (jB.3|) . can be brought into correspondence with the building 
blocks of the non-singlet coefficient function of the polarized structure function gi, A\ 9l ' NS \z) 



displayed in (jB.14J) . For various ER-BL kernels v(x,y), the corresponding DGLAP kernels P(z) 
are obtained by taking the limit ()B.11|) with t = 1 (and i! — > z) taken into account. In Table [7| we 
list some selected results. We mention here that the integration of two "+" forms results again in 
the "+" form, but the contributing terms should be appropriately rearranged. 

Following the procedure explained above, we finally obtain the forward counterparts of the 
elements of the hard-scattering amplitude for the photon-to-pion transition form factor: 



T (0) (cu,x) 



Tp\uj, x) 



5(1 -z), 
-(!-«) + 



3z 1 i / , , 

H ln(l - z) - In , 

1 - z) 1 - z v V ; 



1 + z 



2(1 -z) (1-z) 



31,, . 19 209 (1 + z 2 ) 

— (1 - z) + — ) - y 

12 K J 4(1 -z) 36(1 -z) 



/ 3 / 1 \ 3 
+ 2 (1 " Z) -m^ 



19(l + z 2 )\ H . 
H t r~ hi(l - z) 



(4 

1 + Z 2 



. 5 19(l + z 2 )\ . . , 



1-2 \ 4 

-35(1 - z) . 



5 1 7T' 

+ - | --\n 2 z- -ln 2 (l - z) + 21nzln(l - z) + Li 2 (l - z) + y 



(B.15) 



(B.16) 



(B.17) 



Here, the expressions on the r.h.s. represent the scale independent LO, NLO and /?o-proportional 
NNLO terms of the non-singlet coefficient function (C^ s ) of g\. Similar expressions can be written 
for the terms proportional to ln n (Q 2 / fJ 2 )- Following the notation of [13], the 9(z)9(l—z) factors are 
not shown in Eqs. (|B.15HB.T7j) . We note that the limit of T^(u, 1 — x) (—uj,x) corresponds 

to the antiquark case (C^ s ). 



Our results (jB.15HB.T7j) 7 are numerically in agreement with the Mellin moments and up to a 

typo also with the analytical expression of the n^-proportional term displayed in the appendix of 

Ref. 03]. Namely, in Eq. (A.2) the term 1/3(1 + llz) hzz should read 1/3(1 - lis) In 2. 

7 The representation of the coefficient functions in a form AS(1 — z) + [F(z)] + , as given in Eqs. QB.15IB.f7jl and 
naturally emerging in our calculation, is convenient for the determination of the Mellin moments Cj = z^c(z) 
since the J = Mellin moment corresponds to the term proportional to 5(1 — z). 
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Table 7: Selected forward counterparts of the non-forward quantities. 



6(1 -z) 



(1-3)0(2)0(1-*) 



2z 
1-z 

1 + z 2 

-j— — 0( 
1 — z 

21n(l-z 



z)0(l -z) 
»0(l-z) 
0(z)0(l -*) 



(l-z)(ln2 + 1)0(^)0(1 -z) 



(1-^) + 



1 + z 2 
1-z 



Inz 



0(z)0(l-z) 



IT 



3(1 -z) 



2(1 -z) \nz + ^^\n 2 z 
1 — z 



ln 2 (l-z) n 2Li 2 (l-z) 



1-z 



d(z)d(l-z) 

6(z)6(l-z) 



1-z 



\n(l - u(2y - 1)) 
l-u(2y-l) 



® v(y,x) 



-(1-z) 



2(1 -z) 1-z 



+ ln(l -z) + 



1 + z 2 



1-z 



ln(l — z) — In z 



X0(z)9(l-z) 



49 



C Determination of conformal moments 



In this section we present a method for computing moments with respect to conformal partial 
waves with the index k. We introduce the notation 

(F(x)) k = f dx F(x) C 3 k /2 (2x - 1) , (C.l) 

while iVfc is defined in (J2.1|) . It follows trivially that (F(l — x)) k = (— l) fc (F(x)) k . In the calcu- 
lation of the photon-to-pion transition form factor for the special case \u\ — 1 we encounter the 
functions F(x) of the forms f(x)/x and f(x)/(l — x), with f(x) G {1, ln n (x) ln m (l — x), Li 2 (x), 
Li 3 (x), Si 2 (x) }. 

It is convenient to use the following expression for the Gegenbauer polynomials: 



x(l-x) r 3/2 (2 V] , ]]k 2(2fc + 3) d k fc+1 
The evaluation of the conformal moments, i.e., in our case the evaluation of the expressions 



and 

consists then in calculating the Mellin moments and performing the summation. The Mellin 
moments for the functions we encounter in our calculation are well known (see, for example 
I7d] ) , and most of the nontrivial sums we are left with can be found in jTTj . The sums that 
usually appear are 



n j n 1 

i=l 1 i=l 1 

n n (— 1)* 

1=1 4 i=l 4 



(C.5) 



The functions S m (z) are expressed by ^(z) = dlnT(2;)/dz ones: 

d m 

m = - 7 E + Sxiz - 1) , ^tK*) = [C(m + 1) - S m+1 (z - 1)] . (C.6) 

For (Li 3 (x)/(1 — x)) k and (Si2(x)/(1 — x)) k , the corresponding sums are missing in [ZZj, and to 
obtain them we turn to expressing the relevant functions as convolutions of appropriate functions 



50 



with the known diagonal kernels. Generally, the conformal moments of a kernel [v] + are defined 
by 



vik = ( [v(x,y)]+) 



ik 



J 



dx I dy Cf^ 2 (2x — 1) [v(x, y)] + 



y(i-y) ^3/2 



C?\2y-1), (C.7) 



and for the kernels appearing in this calculation = for I < k and I — k odd. The conformal 
moments of the convolution F(x) = G(y) <g> [v(y,x)] + then take the form 

(F(x))„ = ( G(y) ® [v(y, x)} + ) k = £ (G(y)>, • (C.8) 

«>fc 

As before we use the simplified notation for the diagonal moments 

v k k = v k , (C.9) 

i.e., for the diagonal conformal moments we retain just one index and the relation (|C.8Jl simplifies 
to 

(F (x)) k = ( G(y) ® [v(y, x)} + ) k = (G(y)) k v k . (CIO) 
Hence, in order to determine (Li 3 (x)/(1 — x)) k and (Si 2 (x)/(1 — x)) k , we make use of the identities 



Li 2 (l-y) 
i-y 



® [v a (y,x)\. 



Li 3 (l-x) , C(3) ln(x) ln(x) ln(l - x) 



x 



+ 1 fLi 2 (x) C(2) 



1 - x ' 2(1 -x) 
Li 2 (l-x) C(2)' 



2 \ 1 — x 1 — x 



+ 



x 



(C.ll) 



and 



Li2(y) 



S 12 (x) | C(3) | ln(l-x) | Li 2 (x) 



1 — x 1 — x 
+ C(2) 



x 



2(1 -x) 



ln(l-x) 1 



(C.12) 



1 — x 1 — x J 

The kernels v a and t> are defined in (|3.2j) . (J3.3j) and (|3.9j) . while the corresponding moments can 
be read off from Eqs. (|3.6|) (for e = 0) and (|3.34j) . respectively. 

Finally, in Table |H1 we summarize the conformal moments of the functions relevant to our 
calculation. 

As a byproduct of this calculation, we list the following nontrivial sums: 



E(-i; 



n \ fn + j 

J U + i. 



Sl, 2 tf + 1) 



/I 



E(-i) 



i=0 
rt 



E(-i) J 

i=0 



J fn\ fn + j\ S lt2 (j + 1) 
fn\ fn + j\ S 3 (j + 1) 



j \ j + 1 I j + 1 



n+1 
1 

n + 1 
1 

n + 1 



- (-lWSL 2 (n + 1) + S^ 2 (n - 1)) 
n[n + l) v 7 



■l) n 



n(n + 1) 



S_ 2 (n + l)+£_ 2 (n-l) 



^-^ + (5^ + 1)+^-!)) 



1 



1 



^n(n + l) n 2 (n + l) 2 
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(C.13) 



Table 8: The conformal moments of some relevant functions. 



1 -x 



1 



2(2A; + 3) 
(fc + l)(fc + 2) ~ 2N~ k 



ln(l - 


-x)\ 


1 


1 - 


x Ik 


2N k 


ln 2 (l 


-x)\ 


1 


1 - 


x Ik 


2N k 


ln 3 (l 


-x)\ 


1 


1 - 


x Ik 


2N k 



-Si(fc + 2) -Si(fc) 



Si(A; + 2)-Si(A;) + S 2 (A; + 2) - S 2 (fc) 



Si(fc + 2)-Si(A;)] +2 [_S 3 (fc + 2) -S 3 (fc) 
- Si(fc + 2) - Si(fc)] [S 2 (fc + 2) - S 2 (A;)] } 



+3 



ln(l — x) 



ln 2 (l -x)' 



1 



1 



(k + l)(k + 2) 
-2 



x 



2N k \ (k + l)(k + 2) 



Si(fc + 2)-Si(A;) + l 



ln(x) ln(l — x) 



1 -x 



2iV, 



K(2) - S_ 2 (A; + 2) - S_ 2 (A;) + S 2 (A; + 2) - S 2 (fc) 
- (l - (-1)*) 



+ l)(/s + 2) 



Li 2 (x) 



1 



Li 2 (l-x) 



1 -x 



2iV fc 
1 

2ivT 



C(2) - S 2 (A: + 2) + S 2 (fc) + 



C(2) + S_ 2 (fc + 2) + S_ 2 (fc)- 



(k + l)(k + 2) 
(-1)* 



(fc + l)(fc + 2) 



Lisfx) 



1 — x 



2iV fc 



C(3) 



-1)' 



(A: + 1)(A; + 2) 



C(2) + S_ 2 (A; + 2) + S_ 2 (A;) 



Si 2 (x) 



1 -x 



2N k 



C(3) 



+ 



(A; + !)(*; + 2) 



Si(fc + 2) -Si(fc) 



(fc + l)(fc + 2)- 



Si(A; + 2) - Si(fc)] [S 2 (A; + 2) - S 2 (fc)] } 



52 



which complement the collection of sums found in [77]. 



D Taylor expansions in uj 

We now present the results for the five lowest even partial waves, which are expanded in ui 2 to the 
first seven non-vanishing terms. For brevity, we will not denote the neglected terms. 

The LO result can be simply expanded by means of its representation in terms of hypergeo- 
metrical functions. Employing the identity 



2^1 



'j + l + e,j + 2 + e 



2uj \ = fj/2 + e/2 + 1/2, j/2 + e/2 + 1 
1 + u) 2 j + e + 5/2 



u 2 



(l + uy+ 1 +* 2 'V 2(j + 2 + e) 

and representing the hypergeometrical functions as power series in u, after a few simple manipu- 
lation with T functions we find: 

3 + 2j ~ y/?r(l+j + 2n) /u\ 2n +i 



Tj M -^^ n!r(5/2 +J +n) UJ " (D ^ 
The first few moments read 

T (0) ~ 1 + 0.2c; 2 + 0.0857c; 4 + 0.0476 c; 6 + 0.0303c; 8 + 0.0210c; 10 + 0.0154c; 12 , 
2c; 2 

T 2 (0) ~ — (l + 0.6667c; 2 + 0.4545c; 4 + 0.3263c; 6 + 0.2448c/ + 0.1900c; 10 + 0. 1517c; 12 ) , 
8c; 4 

T 4 (0) ~ fl + 1.1538c; 2 + 1.0769c; 4 + 0.9502c; 6 + 0.8252c; 8 + 0.7152c; 10 + 0.6219c; 12 ) , (D.2) 

315 v ' 

16c; 6 

T 6 (0) ~ — — (l + 1.6471c; 2 + 1.9505c; 4 + 2.0433c; 6 + 2.0211c; 8 + 1.9403c; 10 + 1.8325c; 12 

T (o) 128cj8 f i + 2.1429c; 2 + 3.0745c; 4 + 3.7304c; 6 + 4.1449c; 8 + 4.3736c; 10 + 4.4677c; 
8 109395 V 

The relative error of these approximations is for j = {0, 2, 4, 6, 8} about {0.1%, 0.7%, 2%, 4.4%, 8%} 
for \u\ = 0.8 and increases to {0.6%, 4%, 10%, 19%, 30%} for \u\ = 0.9 . 

The expansion of the Sj(uj) functions from Eq. ()3.63|) can be found in an analogous way: 



= ^n=o^ U,nJ B ,r (5 /2 +J - + n)Uy (D . 3) 

1 ^ 1 v oo r(l+j+2n) (uA 2n ^ ' 

2^ n =0 n!r(5/2+j+n) V 2 / 

<S%» = S 1 (j + 2n) - 51(3/2 +j +n) - S^j) + S 1 (3/2 + j) , (D.4) 
(j, n) = [S® (j, n)] 2 - 5 2 (j + 2n) + S 2 (3/2 + j + n) + S 2 (j) - S 2 (3/2 + j) . 

The approximation of sj(c;) reads 
13c; 2 

4 1} (c;) ~ (l + 0.3642c; 2 + 0.1973c; 4 + 0.1266c; 6 + 0.0894c; 8 + 0.0670c; 10 

53 



oo 2 

s^V) ~ — — (l + 0.3694c/ + 0.2025c/ + 0.1313c/ + 0.0935o> 8 + 0.0707o> 10 ) , (D.5) 
338 ^ ' 

JX)(. ,\ ~ f x + o.3716cj 2 + 0.2048a; 4 + 0.1334c; 6 + 0.0953c/ + 0.0723o> 10 ' 

578 V 



73c; 2 

sj^M ~ fl + 0.3727c; 2 + 0.2059a/ + 0.1344a; 6 + 0.0963c/ + 0.0732a; 1 

294 V 

Note that the prefactor of these series is given by (l — , 5 _^ 2 j)' 2 ) u;2 /4 ~ c//4 and that for j = 
{2, 4, 6, 8} and |o>| = 0.8 the relative error of these approximations is about {1.2%, 2.5%, 4.4%, 6.7%} 
and increases to {6%, 10%, 14%, 19%} for \u\ = 0.9. The approximation of s 2 {ui) is given by 



s {2) < 



sfiu) ~ fl + 14.7002c/ + 10.7004c; 4 + 7.7684c/ + 5.8690c/ + 4.5967c/ 

t \ i 243 ^ 

3c; 2 

sf\u) ~ (l + 44.8125c/ + 33.1180c; 4 + 24.3070c/ + 18.5325c; 8 + 14.6322c/ ) , (D.6) 

2197 ^ 

o 2 

4 2) M ~ f 1 + 100.913c/ + 75.0187a; 4 + 55.3152a; 6 + 42.3458a; 8 + 33.5573a; 1 

6 v ' 4913 v 

c/ 

4 2) (c;) ~ (l + 191.007c/ + 142.411a; 4 + 105.261c/ + 80.7569c/ + 64.1255c/ 

Here we mention that the analytical expansion reads 3c//(5 + 2j) 3 + c/(l + 0(1/(5 + 2 j ) 2 ) / 1 6 . 
Thus, the c/ term is numerically suppressed. In the \u\ — > 1 limit, the functions take the 
values: 

f{uo = 1) = 2^(2j + 3) - S 1 (j + 1) - S 1 (j + 2) - ln(2) , (D.7) 
f (u = 1) ee [sf\u = l)] 2 - AS 2 (2j + 3) + S 2 (j + 2) + S 2 (j + 1) + 2((2) . 

The quantities in the MS scheme are evaluated from Eqs. ()3.16fJ3~2"3*|) . For /i 2 = Q 2 the NLO 
contribution T^\uj , n 2 / Q 2 = 1) reads 

T (1) ~ -2 (l + 0.3333c; 2 + 0.1873c; 4 + 0.1245c; 6 + 0.0904c; 8 + 0.0694c; 10 + 0.0554c; 12 ) , (D.8) 

c; 2 

T 2 {1) ~ — (l + 1.2815c; 2 + 1.0293c; 4 + 0.7770c; 6 + 0.5844c; 8 + 0.4436c; 10 + 0.3407c; 12 ) , 
4184c; 4 

T 4 (1) ~ fl + 1.4076c/ + 1.4571c/ + 1.3683c/ + 1.2370c; 8 + 1.1013a; 10 + 0.9754c/ 2 ) , 

42525 ^ > 

96182c; 6 

T 6 (1) ~ fl + 1.8560c/ + 2.3723c; 4 + 2.6207a> 6 + 2.6949c/ + 2.6644c; 10 + 2.5747a; 12 

2837835 V 

568352c; 8 

T 8 (1) ~ f 1 + 2.3324c; 2 + 3.5542c; 4 + 4.5116a; 6 + 5.1921c/ + 5.6346c; 10 + 5.8894a> ] 

8 57432375 V 

while the conformal moments of the factorization log proportional term, i.e., CpVjTj (u), are 
obtained by multiplying the results (jD.2|) with the values of CpVj, given in Table El At NNLO, 
only the (3q proportional term has been evaluated: 

T$l ~ -3 (l + 0.3505c; 2 + 0.1785c; 4 + 0.1098c/ + 0.0750c/ + 0.0548c; 10 + 0.0420c/ 2 ) , (D.9) 
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T (2) 

PA 



T, 



2) _ 4369c_ / gi59a; 2 + . 7298c 4 + 0.5740a/ + 0.4566c 8 + 0.3691c 10 + 0.3031c 
8640 v 

2356859c 4 / , R s in 10 \ 

(1 + 1.3588c 2 + 1.4063c 4 + 1.3333c 6 + 1.2209c 8 + 1.1022c 10 + 0.9899c 12 ) , 

8505000 v / ' 

(S ) 20352710029c 6 , „ , . „ in 15A 

fe - 222486264000 v 1 + 1M08uJ + 2 -3592c 4 + 2.6233c 6 + 2.7190c 8 + 2.7108c 10 + 2.6415c 12 ) , 

Tf\ ~ 363 f° 06 ° 68 ^ 8 f 1 + 2.3308c 2 + 3.5655a; 4 + 4.5516c 6 + 5.2715c 8 + 5.7582c 10 + 6.0576c 12 ) . 
1 ' 13676945782500 v ' 

To restore the factorization log, one needs 
c 2 

ug ~ — - (l + 0.5657c 2 + 0.3830c 4 + 0.2837c 6 + 0.2220c 8 + 0.1804c 10 ) , (D.10) 
83c 2 

w?2 ~ f 1. + 0.0488c 2 + 0.0308c 4 + 0.0225c 6 + 0.0177c 8 + 0.0146c 10 + 0.0123c 12 ) , 

ft2 216 v ' 

7783c 4 

«L ~ f 1. + 0.0261c 2 + 0.0172c 4 + 0.0130c 6 + 0.0105c 8 + 0.0088c 10 + 0.0076c 12 ) , 

PA 70875 V ' 

~ — fl. + 0.0174c 2 + 0.0118c 4 + 0.0091c 6 + 0.0075c 8 + 0.0064c 10 + 0.0056c 12 ) , 

p ' e 132432300 V / ' 

s _^6991788c_ / Q .0130c 2 + 0.0089c 4 + 0.0070c 6 + 0.0058c 8 + 0.0050c 10 + 0.0044c 12 ) . 
A8 10854718875 V > 

In the CS scheme, the NLO result T° s(1) (c, ^/Q 2 = 1) we find from Eq. (EH)2l : 
T cs(i) _ _ 2 ^ + Q 2uJ 2 + 0.0857c 4 + 0.0476c 6 + 0.0303c 8 + 0.0210c 10 + 0.0154c 12 ) , (D.ll) 

2 

T 2 CS(1) ~ — (l + 1.4691c 2 + 1.2818c 4 + 1.0443c 6 + 0.8466c 8 + 0.6933c 10 + 0.5752c 12 ) , 

cs(i) _ 4184c_ / 1 i 4102a; 2 + ! 4674a; 4 + 1 3875a; 6 + L2 644c 8 + 1.1354c 10 + 1.0148c 12 ) , 
42525 ^ ' 

cs(i) _ 96182c / l g373a; 2 + 2.3344c 4 + 2.5697c 6 + 2.6372c 8 + 2.6050c 10 + 2.5171c 12 ) , 
2837835 ^ ' 

cs(i) _ 568352c A + 2.3052c 2 + 3.4829c 4 + 4.3926c 6 + 5.0299c 8 + 5.4373c 10 + 5.6655c 12 ) . 
8 57432375 V > 

while the factorization log proportional term is the same as in the MS scheme. The NNLO 

correction T^ CS ^(c, [if/Q 2 = l,fi r /Q 2 = 1) for f3 = reads 

T CS(2) ~ 3.6667 (l + 0.2c 2 + 0.0857c 4 + 0.0476c 6 + 0.0303c 8 + 0.021c 10 + 0.0154c 12 ) , (D. 12) 
T cs(2) _ _ 3 2 33kc 2 fl + 0.6769 c 2 + 0.4557c 4 + 0.3224c 6 + 0.2386c 8 + 0.1832c 10 + 0.1448c 12 ) 

1 V / 1/30=0 

T cs(2) _ _o.9338c/ fl + 1.0967c 2 + 0.9764c 4 + 0.8272c 6 + 0.6939c 8 + 0.5838c 10 + 0.4949c 12 ) 
4 V /|A)=o 

T <cs(2) ~ -0.21c 6 fl + 1.5133c 2 + 1.6617c 4 + 1.6285c 6 + 1.5184c 8 + 1.3829c 10 + 1.2459c 12 ) 

T cs(2) _ _ .0427c 8 fl + 1.905c 2 + 2.4505c 4 + 2.6858c 6 + 2.7134c 8 + 2.6181c 10 + 2.4577c 12 ) 

8 v / |A)=o 

The Q 2 independent and /^-proportional term in the CS scheme is the same as in the MS one, 
given in Eq. (|D.12|) . The factorization and renormalization log proportional terms, appearing in 
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Table 9: The first four non-vanishing and even eigenvalues of the evolution kernel ()1.5|) to NNLO 

_____ /q\ 

accuracy (for color decomposition see ([3.39)) ) . The values of Vj are taken from the non-singlet 
result for the deep inelastic structure function F 3 |5H] . 



3 


CpVj 


/-f2 Up 

C F v F , + 2N v ad 




vf = -l/2 7 f 


2 


-2. mis 


3.41307 


-2.88194 


-155.614 + 24.5592/1/ + 0.220250 n) 


4 


-A.OAAAA 


7.15867 


-4.32389 


-215.118 + 34.7698n/ + 0.295776n| 


6 


-4.89048 


9.82554 


-5.30857 


-254.562 + 41.3602 n f + 0.342420 n) 


8 


-5.52910 


11.86905 


-6.06196 


-284.650 + 46.3238 n f + 0.375806 n) 



Eq. (|3.64|) can be easily restored by means of the results from Table |H1 Eqs. ()D.2ID.5|) as well as 
Eqs. lj3ireilD.8ID.10ID.llll . 

The difference between the CS and CS schemes arises only from the /3o-proportional terms 
(J3.69j) . The /3o-proportional NNLO term Tpj (u) reads 

T^o (2) ~ -3 (l + 0.2c; 2 + 0.0857c/ + 0.0476c; 6 + 0.0303c; 8 + 0.0210c; 10 + 0.0154c; 12 ) , (D.13) 
Tpf 2) ~ 0.50567c; 2 (l. + 0.8496c; 2 + 0.6431c; 4 + 0.4900c; 6 + 0.3820c; 8 + 0.3048c; 10 + 0.2482c; 12 ) , 
Tpf 2) ~ 0.277114c; 4 (l + 1.2512c; 2 + 1.2252c; 4 + 1.1162c; 6 + 0.9919c; 8 + 0.8747c; 10 + 0.7710c; 12 ) , 
Tpf 2) ~ 0.0914785c; 6 (l. + 1.7236c; 2 + 2.1049c; 4 + 2.2550c; 6 + 2.2689c; 8 + 2.2076c; 10 + 2.1078c; 12 ) , 
Tpf 2) ~ 0.02656c; 8 (l + 2.2092c; 2 + 3.2413c; 4 + 4.0009c; 6 + 4.5064, c; 8 + 4.8081c; 10 + 4.9570c; 12 ) . 

The restoration of the factorization log in the (3q sector requires the knowledge of vpj, given in 
Table M 

E Reconstruction of the hard-scattering amplitude in the 
momentum fraction representation 

Let us now discuss the reconstruction of the hard-scattering amplitude in the momentum fraction 
representation from the conformal moments. This technical problem is of immense importance for 
the discussion of two-photon processes in which the operator product expansion is not convergent. 
The solution is known in forward kinematics and is given by the Mellin transformation of moments. 
In non-forward kinematics, the problem is solved in principle however, one has to evaluate 
rather cumbersome integrals. Here we propose a simple approximative solution which is based on 
the asymptotic behaviour, presented in Section 4.1, and it is applicable to deeply virtual Compton 
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scattering. In the following we demonstrate this method at NLO order and its generalization to 
higher orders is straightforward. 

First, we consider the quality of the asymptotic formulas (j4.1j) . Since the neglected terms are 
of the order + 2), we expect that the approximations (|4.1|) have an accuracy of the level of 
10% for j = 10. Surprisingly, the accuracy is already below 1% in both cases, which indicates that 
the + 2) term is small. Thus, we completely neglect such terms and improve the approximate 
formula by adding subleading terms of the form 

(j + 1)(j + 2) + ( J + l)2( J + 2)2 + -"' 

where the coefficients a, ... ,5 are determined from a fit of the lowest moments. In this way, 
we obtain an approximation that is better than 1% for all moments. Now we can recover an 
approximate expression for the hard-scattering amplitude as a convolution by the following recipe: 

• Substitute the LO Wilson-coefficients by the corresponding hard-scattering amplitude: 

2j + 3 ) 1 
(j + l)(j + 2) ~* 2(1-*)' 

• Restore kernels: 

cons cons I, 1 -> v a (x,y), Si(l + j) -> --[v b (x, y)}+ + 1 1. 

U + lJU + 2) 2 

• Consider the multiplication of the above given conformal moments, which corresponds to 
convolution in the momentum fraction space. 

Here we have introduced a shorthand notation for the identity I = 5(x — y). In the MS scheme, 

using Eq. ()3.32|) and the recipe given above we restore the exact expression ()3.28|) in the momentum 

fraction space. For the CS scheme, using the improved form (|4.1|l . we get a good approximation 

of Eq. ()3.45|) . for the lowest moments also, by taking 

83 73 
a = — + 8C(2) - 26 ln(2) « 3.438, 7 = -— - 12((2) + 48 ln(2) « -1.068, (3 = 5 = (E.l) 

-LU 

and the hard-scattering part reads 

T£ s{1) (x) « (2N+- 6.273 1) — 2.952 1 + 3.438 v a — 1.068 v a <8> u a | (x) 

_ 0.25 ln 2 (l - x) - 2.136 ln(l - x) - 8.224 ln(l - x) [0.25 ln(l - x) - 6.642] 
~ 2(1 -x) + Yx 

+ 1.068 Li 'g 1 -_ L ^ 1 » . (E.2) 

After analytical continuation in x this result corresponds to the NLO correction of the deeply 
virtual Compton scattering for the quark-quark channel in the parity odd sector. 
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